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Abstract. We present a theory and applications of discrete exterior calculus on simplicial complexes of 
arbitrary finite dimension. This can be thought of as calculus on a discrete space. Our theory includes not 
only discrete differential forms but also discrete vector fields and the operators acting on these objects. This 
allows us to address the various interactions between forms and vector fields (such as Lie derivatives) which 
are important in applications. Previous attempts at discrete exterior calculus have addressed only differential 
forms. We also introduce the notion of a circumcentric dual of a simplicial complex. The importance of dual 
complexes in this field has been well understood, but previous researchers have used barycentric subdivision 
or barycentric duals. We show that the use of circumcentric duals is crucial in arriving at a theory of discrete 
exterior calculus that admits both vector fields and forms. 
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1. Introduction 

This work presents a theory of discrete exterior calculus (DEC) motivated by potential applications 
in computational methods for field theories such as elasticity, fluids, and electromagnetism. In addition, it 
provides much needed mathematical machinery to enable a systematic development of numerical schemes 
that mirror the approach of geometric mechanics. 

l 
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This theory has a long history that we shall outline below in but we aim at a comprehensive, systematic, 
as well as useful, treatment. Many previous works, as we shall review, are incomplete both in terms of the 
objects that they treat as well as the types of meshes that they allow. 

Our vision of this theory is that it should proceed ab initio as a discrete theory that parallels the c ontinuous 
one. General views of the subject area of DEC are common in the literature (see, for instance, iMattiussl 
2000J), but they usually stress the process of discretizing a continuous theory and the overall approach is 
tied to this goal. However, if one takes the point of view that the discrete theory can, and indeed should, 
stand in its own right, then the range of application areas naturally is enriched and increases. 

Convergence and consistency considerations alone are inadequate to discriminate between the various 
choices of discretization available to the numerical analyst, and only by requiring, when appropriate, that 
the discretization exhibits discrete analogues of continuous properties of interest can we begin to address the 
question of what makes a discrete theory a canonical discretization of a continuous one. 

Applications to Variational Problems. One of the major application areas we envision is to variational 
problems, be they in mechanics or optimal control. One of the key ingredients in this direction that we 
imagine will play a key role in the future is that of AVI' s (asynchronous v ariational integrators) designed for 
the numerical integration of mechanical systems, as in iLew et al . 2003]. These are integration algorithms 
that respect some of the key features of the continuous theory, such as their multi-symplectic nature and 
exact conservation laws. They do so by discretizing the underlying variational principles of mechanics rather 
than discretizing the equations. It is well-known (see the reference just mentioned for some of the literature) 
that variational problems come equipped with a rich exterior calculus structure and so on the discrete level, 
such structures will be enhanced by the availability of a discrete exterior calculus. One of the objectives of 
this chapter is to fill this gap. 

Structured Constraints. There are many constraints in numerical algorithms that naturally involve 
differential forms, such as the divergence constraint for incompressibility of fluids, as well as the fact that 
differential forms are naturally the fields in electromagnetism, and some of Maxwell's equations are expressed 
in terms of the divergence and curl operations on these fields. Preserving, as in the mimetic differencing 
literature, such features directly on the discrete level is another one of the goals, overlapping with our goals 
for variational problems. 

Lattice Theories. Periodic crystalline lattices are of important practical interest in material science, and the 
anisotropic nature of the material properties arises from the geometry and connectivity of the intcrmolccular 
bonds in the lattice. It is natural to model these lattices as inherently discrete objects, and an understanding 
of discrete curvature that arises from DEC is particularly relevant, since part of the potential energy arises 
from stretched bonds that can be associated with discrete curvature in the underlying relaxed configuration of 
the lattice. In particular, this could yield a more detailed geometric understanding of what happens at grain 
boundaries. Lattice defects can also be associated with discrete curvature when appropriately interpreted. 
The introduction of a discrete notion of curvature will lay the foundations for a better understanding of the 
role of geometry in the material properties of solids. 

Some of the Key Theoretical Accomplishments. Our development of discrete exterior calculus in- 
cludes discrete differential forms, the Hodge star operator, the wedge product, the exterior derivative, as 
well as contraction and the Lie derivative. For example, this approach leads to the proper definition of 
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discrete divergence and curl operators and has already resulted in applica tions like a discrete Hodge type 
decomposition of 3D vector fields on irregular grids — see iTong et al.l ,20031. 

Context. We present the theory and some applications of DEC in the context of simplicial complexes of 
arbitrary finite dimension. 

Methodology. We believe that the correct way to proceed with this program is to develop, as we have 
already stressed, ab initio, a calculus on discrete manifolds which paralle ls the calculus on smooth manifolds 
of arbitrary finite dimension. Chapters 6 and 7 of I Abraham et all 198Sj are a good source for the concepts 
and definitions in the smooth case. However we have tried to make this chapter as self-contained as possible. 
Indeed, one advantage of developing a calculus on discrete manifolds, as we do here, is pedagogical. By 
using concrete examples of discrete two- and three-dimensional spaces one can explain most of calculus on 
manifolds at least formally as we will do using the examples in this chapter. The machinery of Ricmannian 
manifolds and general manifold theory from the smooth case is, strictly speaking, not required in the discrete 
world. The technical terms that are used in this introduction will be defined in subsequent sections, but 
they should be already familiar to someone who knows the usual exterior calculus on smooth manifolds. 

The Objects in DEC. To develop a discrete theory, one must define discrete differential forms along with 
vector fields and operators involving these. Once discrete forms and vector fields are defined, a calculus 
can be developed by defining the discrete exterior derivative (d), codifferential (6) and Hodge star (*) for 
operating on forms, discrete wedge product (A) for combining forms, discrete flat (b) and sharp (jj) operators 
for going between vector fields and 1-forms and discrete contraction operator (ix) for combining forms and 
vector fields. Once these are done, one can then define other useful operators. For example, a discrete Lie 
derivative {£x) can be defined by requiring that the Cartan magic (or homotopy) formula hold. A discrete 
divergence in any dimension can be defined. A discrete Laplace-deRham operator (A) can be defined using 
the usual definition of d£ + £d. When applied to functions, this is the same as the discrete Laplace-Beltrami 
operator (V 2 ), which is the defined as divocurl. We define all these operators in this chapter. 

The discrete manifolds we work with are simplicial complexes. We will recall the standard formal def- 
initions in 321 but familiar examples of simplicial complexes are meshes of triangles embedded in R 3 and 
meshes made up of tetrahedra occupying a portion of R 3 . We will assume that the angles and lengths on 
such discrete manifolds are computed in the embedding space M. N using the standard metric of that space. 
In other words, in this chapter we do not address the issue of how to discretize a given smooth Riemannian 
manifold, and how to embed it in R , since there may be many ways to do this. For example, SO (3) can 
be embedded in K 9 with a constraint, or as the unit quaternions in R 4 . Another potentially important con- 
sideration in discrctizing the manifold is that the topology of the simplicial complex should be the same as 
the mani fold to be discret i zed. T his can be verified using the methods of co mputational homology ( see, fo r 
example. iKaczvnski et all 2004j ). or discrete Morse theory (see, for example, iFormanl |2002j | , IWoodl |2003| ) . 
For the purposes of discrete exterior calculus, only local metric information is required, and we will comment 
towards the end of SjSI how to address the issue of embedding in a local fashion, as well as the criterion for a 
good global embedding. 

Our development in this chapter is for the most part formal in that we choose appropriate geometric 
definitions of the various objects and quantities involved. For the most part, we do not prove that these 
definitions converge to the smooth counterparts. The definitions are chosen so as to make some important 
theorems like the generalized Stokes' theorem true by definition. Moreover, in the cases where previous 
results are available, we have checked that the operators we obtain match the ones obtained by other means, 
such as variational derivations. 
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2. History and Previous Work 

The use of simplicial chains and cochains as th e basic bu i lding blocks for a di s crete exterior ca lculus 
has appeared in several papers. See, for instance, ISen et ail 2000j . lAdamd 1996| . IBossavitl |2002d| . and 
references therein. These a uthors view form s as linearly interpolated versions of smooth differential forms, a 
viewpoint originating from IWhitnevI 1957t | . who introduced the Whitney and deRham maps that establish 
an isomorphism between simplicial cochains and Lipschitz differential forms. 

We will, however, view discrete forms as real-valued linear functions on the space of chains. These are 
inherently discrete objects that can be paired with chains of oriented simplices, or their geometric duals, 
by the bilinear pairing of evaluation. In the next chapter, where we consider applications involving the 
curvature of a discrete space, we will relax the condition that discrete forms are real-valued, and consider 
group- valued forms. 

Intuitively, this natural pairing of evalu ation can be tho ught of as integra tion of the discrete form over 
the chain. This difference from the work of lSen et alJ |200n| and lAdamsl |l99fi| is apparent in the definitions 
of operations like the wedge product as well. 

There is also much interest in a discrete exte ri or calculus in the computa t ional electromagnetism commu - 



1 here is also much interest m a discrete exterior calculus m the computational electromagnetism commu- 
nity, as represented by IBossavitl l200lil2002albO.lGross and Kotiugal l2nf)ll.lHiptmairl |l999ll2001alhl l2002j | , 
iMattiussil |l997ll20odj . lNicolaides and Wanel ll99Sl. iTeixeiral l2Q0ll. andlTontil l2002j . 

Many of the authors cited above, for example. IBossavitl |20fl2cj . lSen et al.l |20f)f)j . and lHiotmairl |20fl2l | , also 
introduce the notions of dual complexes in order to construct the Hodge star operator. With the exception 
of Hiptmair, they use barycentric duals. This works if one develops a theory of discrete forms and does not 
introduce discrete vector fields. We show later that to introduce discrete vector fields into the theory the 

notion of circumcentric duals seems to be important. 

Other authors, such as lMoritd [200Cj . iMoritz and Schwalml |20f)lj . ISchwalm et alJ [l999j , have incorporated 



vector fields into the cochain based approach to exterior calculus by identifying vector fields with cochains, 
and having them supported on the same mesh. This is ultimately an unsatisfactory approach, since dual 
meshes are essential as a means of encoding physically relevant phenomena such as fluxes across boundaries. 

The use of primal and dual meshes arises most often as staggered meshes in finite volume and finite 
difference methods. In fluid computations, for example, the density is often a cell-centered quantity, which 
can either be represented as a primal object by being associated with the 3-cell, or as a dual object associated 
with the 0-cell at the center of the 3-cell. Similarly, the flux across boundaries can be associated with the 
2-cells that make up the boundary, or the I-cell which is normal to the bound ary. 

Another approach to a discrete exterior calculus is presented in lDezinl |1995j . He defines a one-dimensional 
discretization of the real line in much the same way we would. However, to generalize to higher dimensions 
he introduces a tensor product of this space. This results in logically rectangular meshes. Our ca lculus, 
howev er, is defined over simplicial meshes. A further difference is that like other authors in this field, iDezinl 
1995] does not introduce vector fields into his theory. 

A r elated effort for three-dimensional domains with logically rectangular meshes is that o flMansfield and HvdonI 
who established a variational complex for difference equations by constructing a discrete homotopy op- 
erator. We construct an analogous homotopy operator for simplicial meshes in proving the discrete Poincare 
lemma. 



3. Primal Simplicial Complex and Dual Cell Complex 

In constructing the discretization of a continuous problem in the context of our formulation of discrete 
exterior calculus, we first discretize the manifold of interest as a simplicial complex. While this is typically 
in the form of a simplicial complex that is embedded into Euclidean space, it is only necessary to have an 
abstract simplicial complex, along with a local metric defined on adjacent vertices. This abstract setting will 
be addressed further toward the end of this section. 
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We will now recall some basic definitions of simplices and simplicial co mplexes, which a re standard from 
simplicial algebraic topology. A more extensive treatment can be found in 

Definition 3.1. A k-simplex is the convex span of k + 1 geometrically independent points, 

( k 



a k = [vo, vi,...,Vk] = \ a l Vi 



a 1 > 0, a * 

i=0 



k i=0 

The points VQ,...,Vk are called the vertices of the simplex, and the number k is called the dimension of 
the simplex. Any simplex spanned by a (proper) subset of {vo, . . . ,Vk} is called a (proper) face of a k . If 
a 1 is a proper face of a k , we denote this by a 1 -< a k . 

Example 3.1. Consider 3 non-collinear points i>o, «i and V2 in M 3 . Then, these three points individually 
are examples of 0-simplices, to which an orientation is assigned through the choice of a sign. Examples of 
1-simplices are the oriented line segments [vo,v±], [vi,V2] and [uojVa]. By writing the vertices in that order 
we have given orientations to these 1-simplices, i.e., [vq,Vi] is oriented from vq to v\. The triangle [vq, v±, V2] 
is a 2-simplex oriented in counterclockwise direction. Note that the orientation of [vq,V2] does not agree with 
that of the triangle. 

Definition 3.2. A simplicial complex K in R N is a collection of simplices in M. N , such that, 

(1) Every face of a simplex of K is in K . 

(2) The intersection of any two simplices of K is a face of each of them. 

Definition 3.3. A simplicial triangulation of a polytope \K\ is a simplicial complex K such that the 
union of the simplices of K recovers the polytope \K\. 

Definition 3.4. If L is a subcollection of K that contains all faces of its elements, then L is a simplicial 
complex in its own right, and it is called a subcomplex of K . One subcomplex of K is the collection of all 
simplices of K of dimension at most k, which is called the k-skeleton of K, and is denoted K^ k \ 

Circumcentric Subdivision. We will also use the notion of a circumcentric dual or Voronoi mesh of the 
given primal mesh. We will point to the importance of this choice later on in |Qand|ni We call the Voronoi 
dual a circumcentric dual since the dual of a simplex is its circumcenter (equidistant from all vertices of the 
simplex) . 

Definition 3.5. The circumcenter of a k-simplex a k is given by the center of the k-circumsphere, where 
the k-circumsphere is the unique k-sphere that has all k + 1 vertices of a k on its surface. Equivalently, 
the circumcenter is the unique point in the k-dimensional affine space that contains the k-simplex that is 
equidistant from all the k + 1 nodes of the simplex. We will denote the circumcenter of a simplex a k by 
c(a k ). 

The circumcenter of a simplex o~ k can be obtained by taking the intersection of the normals to the (k — 1)- 
dimensional faces of the simplex, where the normals are emanating from the circumcenter of the face. This 
allows us to recursively compute the circumcenter. 

If we are given the nodes which describe the primal mesh, we can construct a simplicial triangulation by 
using the Delaunay triangulation, since this ensures that the circumcenter of a simplex is always a point 
within the simplex. Otherwise we assume that a nice mesh has been given to us, i.e., it is such that the 
circumcenters lie within the simplices. While this is not be essential for our theory it makes some proofs 
simpler. For some computations the Delaunay triangulation is desirable in that it reduces the maximum 
aspect ratio of the mesh, which is a factor in determining the rate at which the corresponding numerical 
scheme converges. Bu t in practice there are many problems for which Delaunay triangulations are a bad 
idea. See, for example. ISchewchuckl 2002) . We will address such computational issues in a separate work. 



6 



MATHIEU DESBRUN, ANIL N. HIRANI, MELVIN LEOK, AND JERROLD E. MARSDEN 



Definition 3.6. The circumcentric subdivision of a simplicial complex is given by the collection of all 
simplices of the form 

[c(cr ), . . . ,c(cr fe )], 

where er -< a\ ~< . . . -< a^, or equivalently, that o~i is a proper face of o~j for all i < j . 

Circumcentric Dual. Wc construct a circumcentric dual to a fc-simplex using the circumcentric duality 
operator, which is introduced below. 

Definition 3.7. The circumcentric duality operator is given by 

*(<7 fc )= Yl e ak _ an [c{a k ),c{a k+l ),...,c(a n )}, 

a k -<a k + 1 ^ l ...-<<j n 

where the e a k a n coefficient ensures that the orientation of [c(a k ), c(a k+1 ),..., c(cr™)l is consistent with 
the orientation of the primal simplex, and the ambient volume-form. 

Orienting a k is equivalent to choosing a ordered basis, which we shall denote by dx 1 A ... A dx k . Similarly, 
[c(<7 fc ), c(f7 fc+1 ), . . . , c(er™)] has an orientation denoted by dx k+1 A ... A dx n . If the orientation corresponding 
to dx 1 A ... A dx n is consistent with the volume-form on the manifold, then e a k = 1, otherwise it takes 
the value —1. 

We immediately see from the construction of the circumcentric duality operator that the dual elements 
can be realized as a submesh of the first circumcentric subdivision, since it consists of elements of the form 
[c((Jo), • • ■ , c(cr/c)], which are, by definition, part of the first circumcentric subdivision. 

Example 3.2. The circumcentric duality operator maps a 0-simplex into the convex hull generated by the 
circumcenters of n- simplices that contain the 0-simplex, 



,(<7°) = {5> ffn cK) 



a a n > 0,^a ff ™ = 1 



and the circumcentric duality operator maps a n— simplex into the circumcenter of the n— simplex, 

*(0 = c(a"). 

This is more clearly illustrated in Figure^ where the primal and dual elements are color coded to represent 
the dual relationship between the elements in the primal and dual mesh. 






(a) Primal 



(b) Dual 



(c) First subdivision 



Figure 1. Primal, and dual meshes, as chains in the first circumcentric subdivision. 



The choice of a circumcentric dual is significant, since it allows us to recover geometrically important 
objects such as normals to (n — l)-dimensional faces, which are obtained by taking their circumcentric dual, 
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whereas, if we were to use a barycentric dual, the dual to a (n — l)-dimensional face would not be normal to 
it. 



which is represented by 
• A (vk — Ufe-i), which we 



Orientation of the Dual Cell. Notice that given an oriented simplex <j k , 
[vq, ...,Vk], the orientation is equivalently represented by (v% — vq) A (i>2 — v%) A . 
denote by, 

[vo, ■ ■ ■ ,Vk] ~ (vi - v ) A (v 2 - Vi) A . . . A (Vk ~ Ufc-l), 
which is an equivalence at the level of orientation. It would be nice to express our criterion for determining 
the orientation of the dual cell in terms of the (fc + l)-vertex representation. 

To determine the orientation of the (n — fc)-simplex given by [c(a k ), c(a k+1 ), . . . , c(a n )], or equivalently, 
dx k+1 A ... A dx n , we consider the n-simplex given by [c(cr°), . . . , c(a n )], where a ~< . . . -< a k . This is related 
to the expression dx 1 A ... A dx n , up to a sign determined by the relative orientation of [c(er ), . . . , c(a k )] 
and a k . Thus, we have that 

dx 1 A ... A dx n ~ sgn([c(a°), . . . , c(a k )}, a k )[c(a°), c(a n )} . 

Then, we need to check that dx 1 A. . .Adx n is consistent with the volume-form on the manifold, which is repre- 
sented by the orientation of a n . Thus, we have that the correct orientation for the [c(a k ), c(a k+1 ), . . . , c(cr™)] 
term is given by, 

sgn([c(a°), . . . , c(a% a k ) ■ Sgn([c(a°), . . . , c{a n )],a n ). 
These two representations of the choice of orientation for the dual cells are equivalent, but the combinatorial 
definition above might be preferable for the purposes of implementation. 

Example 3.3. We would like to compute the orientation of the dual of a 1- simplex, in two dimensions, 
given the orientation of the two neighboring 2-simplices. 

Given a simplicial complex, as shown in Figure 2(a)\ we consider a 2-simplex of the form [c(cr°), c(ct 1 ), c(er 2 )], 



which is illustrated in Figure 2(b) 









(a) Simplicial complex 



(b) 2-simplex 



(c) ^cr 1 



Figure 2. Orienting the dual of a cell. 



Notice that the orientation is consistent with the given orientation of the 2-simplex, but it is not consistent 
with the orientation of the primal 1-simplex, so the orientation should be reversed, to give the dual cell 



illustrated in Figure 2(c) 



We summarize the results for the induced orientation of dual cells for the other 2-simplices of the form 
^(cr ),^ 1 ),^ 2 )], in Tabled 



Orientation of the Dual of a Dual Cell. While the circumcentric duality operator is a map from the 
primal simplicial complex to the dual cell complex, we can formally extend the circumcentric duality operator 
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Table 1. Determining the induced orientation of a dual cell. 



[c(a%c(a'),c(a 2 )] 




4\ 






V 


\d/ 


V 




S gn([c(a°),c(a%a l ) 




+ 


+ 




sgnCKa ),^ 1 ),^ 2 )],^) 


+ 




+ 




sgn([c(a ) >C (a%*i) 
•sgnCKa ),^ 1 ),^ 2 )],^) 










V 


V 







to a map from the dual cell complex to the primal simplicial complex. However, we need to be slightly careful 
about the orientation of primal simplex we recover from applying the circumcentric duality operator twice. 

We have that, ** (a k ) — ±er fe , where the sign is chosen to ensure the appropriate choice of orientation. 
If, as before, a k has an orientation represented by dx 1 A ... A dx k , and *<r fc has an orientation represented 
by dx k+1 A ... A dx n , then the orientation of ** (a k ) is chosen so that dx k+1 A ... A dx 11 A da; 1 A ... A dx k is 
consistent with the ambient volume-form. Since, by construction, *(<r fe ), dx 1 A ... A dx n has an orientation 
consistent with the ambient volume-form, we need only compare dx k+1 A ... A dx n A dx 1 A ... A dx k with 
dx 1 A ... A dx n . Notice that it takes n — k transpositions to get the dx 1 term in front of the dx k+1 A ... A dx n 
terms, and we need to do this k times for each term of dx 1 A ... A dx k , so it follows that the sign is simply 
given by (— l) fe (™~ fc ) ) or equivalently, 

(3.1) ★★(cr fc ) = (-l) fc ("- fc V fc . 

A similar relationship holds if we use a dual cell instead of the primal simplex o~ k . 

Support Volume of a Primal Simplex and Its Dual Cell. We can think of a cochain as being 
constructed out of a basis consisting of cosimplices or cocells with value 1 on a single simplex or cell, and 
otherwise. The way to visualize this cosimplex is that it is associated with a differential form that has 
support on what we will refer to as the support volume associated with a given simplex or cell. 

Definition 3.8. The support volume of a simplex o~ k is a n-volume given by the convex hull of the geometric 
union of the simplex and its circumcentric dual. This is given by 

V a k = convexhull(cr fe ,*cr fe ) n \K\. 

The intersection with \K\ is necessary to ensure that the support volume does not extend beyond the 
polytope \K\ which would otherwise occur if \K\ is nonconvex. 

We extend the notion of a support volume to a dual cell ★tr fc by similarly defining 

y*a k = convexhull(*cr' £ ,**o' fe ) PI \K\ = V a k. 
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To clarify this definition, we will consider some examples of simplices, their dual cells, and their corre- 
sponding support volumes. For two-dimensional simplicial complexes, this is illustrated in Table 

Table 2. Primal simplices, dual cells, and support volumes in two dimensions. 



Primal Simplex 


Dual Cell 


Support Volume 




A 


m 


t 


(T°, 0-simplex 


*cr°, 2-cell 


V a o = 












er 1 , 1-simplex 


★cr 1 , 1-cell 


V a x = 






• 


A 


k 


er 2 , 2-simplex 


*<t 2 , 0-cell 


V* = 





The support volume has the nice property that at each dimension, it partitions the polytope \K\ into 
distinct non-intersecting regions associated with each individual fc-simplex. For any two distinct fc-simplices, 
the intersection of their corresponding support volumes have measure zero, and the union of the support 
volumes of all fc-simplices recovers the original polytope \K\. 

Notice, from our construction, that the support volume of a simplex and its dual cell are the same, which 
suggests that there is an identification between cochains on fc-simplices and cochains on (n — fc)-cells. This 
is indeed the case, and is a concept associated with the Hodge star for differential forms. 

Examples of simplices, their dual cells, and the corresponding support volumes in three dimensions are 
given in Table [3] 

In our subsequent discussion, we will assume that we are given a simplicial complex K of dimension n in 
R . Thus, the highest-dimensional simplex in the complex is of dimension n and each 0-simplex (vertex) is 
in R N . One can obtain this, for example, by starting from O-simplices, i.e., vertices, and then constructing a 
Delaunay triangulation, using the vertices as sites. Often, our examples will be for two-dimensional discrete 
surfaces in R 3 made up of triangles (here n = 2 and N = 3) or three-dimensional manifolds made of 
tetrahedra, possibly embedded in a higher-dimensional space. 

Cell Complexes. The circumcentric dual of a primal simplicial complex is an example of a cell complex. 
The definition of a cell complex follows. 

Definition 3.9. A cell complex *K in R N is a collection of cells in 1& N such that, 

(1) There is a partial ordering of cells in *K, a k -< a , which is read as & k is a face of a 1 . 

(2) The intersection of any two cells in -kK , is either a face of each of them, or it is empty. 

(3) The boundary of a cell is expressible as a sum of its proper faces. 
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Table 3. Primal simplices, dual cells, and support volumes in three dimensions. 



Primal Simplex 



Dual Cell 



Support Volume 




a , 0-simplex 




*cr°, 3-cell 




cr 1 , 1-simplex 




*<t\ 2-cell 





cr 2 , 2-simplex 




*(7 2 1-cell 




V a 2 = 






cr 3 , 3-simplex *er 3 , 0-cell 



We will see in the next section that the notion of boundary in the circumcentric dual has to be modified 
slightly from the geometric notion of a boundary in order for the circumcentric dual to be made into a cell 
complex. 



4. Local and Global Embeddings 

While it is computationally more convenient to have a global embedding of the simplicial complex into a 
higher-dimensional ambient space to account for non-flat manifolds it suffices to have an abstract simplicial 
complex along with a local metric on vertices. The metric is local in the sense that distances between two 
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vertices are only denned if they are part of a common rt-simplex in the abstract simplicial complex. Then, 
the local metric is a map d : {(vo, «i) | vo,vi G K^°\ [vo,vi] -< er™ G K} — > M. 
The axioms for a local metric are as follows, 

Positive.: d(v ,vi) > 0, and d(vo,v n ) = 0, V[uo,i;i] -< cr n G K. 

Strictly Positive.: If d(v ,vi) ~ 0, then v — vi, y[v ,vi] -< cr™ G -ftT. 

Symmetry.: d(v ,vi) = d(vi,v a ), V[u , «i] -< a n G K. 

Triangle Inequality.: d(v ,v 2 ) < d(v Q ,vi) + d(vi,v 2 ), V[vq,Vi,V2] -< cr n G K. 
This allows us to embed each n-simplex locally into R ra , and thereby compute all the necessary metric 
dependent quantities in our formulation. For example, the volume of a A;-dual cell will be computed as the 
sum of the k- volumes of the dual cell restricted to each n-simplex in its local embedding into W 1 . 

This notion of local metrics and local embeddings is consistent with the point of view that exterior calculus 
is a local theory with operators that operate on objects in the tangent and cotangent space of a fixed point. 
The issue of comparing objects in differen t tangent spaces is addressed in the discrete theory of connections 



This also provides us with a criterion for evaluating a global embedding. The embedding should be such 
that the metric of the ambient space M. N restricted to the vertices of the complex, thought of as points in WL N , 
agrees with the local metric imposed on the abstract simplicial complex. A global embedding that satisfies 
this condition will produce the same numerical results in discrete exterior calculus as that obtained using 
the local embedding method. 

It is essential that the metric condition we impose is local, since the notion of distances between points in 
a manifold which are far away is not a well-defined concept, nor is it particularly useful for embeddings. As 
the simple example below illustrates, there may not exist any global embeddings into Euclidean space that 
satisfies a metric constraint imposed for all possible pairs of vertices. 

Example 4.1. Consider a circle, with the distance between two points given by the minimal arc length. 
Consider a discretization given by 4 equidistant points on the circle, labelled vq, . . . ,v^, with the metric 
distances as follows, 



where the indices are evaluated modulo 4? an d this distance function is extended to a metric on all pairs of 
vertices by symmetry. It is easy to verify that this distance function is indeed a metric on vertices. 




d(v l7 v i+1 ) = l,d(vi,Vi +2 ) = 2, 




® 



If we only use the local metric constraint, then we only require that adjacent vertices are separated by 1, and 
the following is an embedding of the simplicial complex into R 2 , 
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// ; however, we use the metric defined on all possible pairs of vertices, by considering vq, v i, V2, we have that 
d(vo,vi) + d{v\,V2) = d(vo,V2). Since we are embedding these points into a Euclidean space, it follows that 
Vq,vi,V2 are collinear. 

Similarly, by considering 1)0,1)2,1)3, we conclude that they are collinear as well, and that 1)1,1)3 are coinci- 
dent, which contradicts d(i)i, 1)3) = 2. Thus, we find that there does not exist a global embedding of the circle 
into Euclidean space if we reguire that the embedding is consistent with the metric on vertices defined for all 
possible pairs of vertices. 



5. Differential Forms and Exterior Derivative 

We will now define discrete differential forms. We will use some terms (which we will define) from algebraic 
topology, but it will become clear by looking at the examples that one can gain a clear and working notion 
of what a discrete form is without any al gebraic topology . We start with a few definitions for which more 
details can be found on page 26 and 27 of iMunkred |l984 |- 

Definition 5.1. Let K be a simplicial complex. We denote the free abelian group generated by a basis 
consisting of oriented k-simplices by Ck (K;"L) . This is the space of finite formal sums of the k-simplices, 
with coefficients inl. Elements of Ck(K;Z) are called k-chains. 

Example 5.1. Figure^ shows examples of 1-chains and 2 -chains. 




1— chain 2— chain 



Figure 3. Examples of chains. 



We view discrete /c-forms as maps from the space of fc-chains to R. Recalling that the space of £;-chains 
is a group, we require that these maps be homomorphisms into the additive group R. Thus, discrete forms 
are what are called cochains in algebraic topology. We will define cochains below in the definition of forms 
but fo r more context and more details readers can refer to any algebraic topology text, for example, page 
251 of lMunkresI |l98^ . 

This point of view of f orms as cocha i ns is not new . The idea of de fining forms as c ochai ns appears, for 
example, in the works of lAdamd |l996f . iDezin] |l995j . iHiptmairl |l999| . and lSen et ail |200tjj . Our point of 
departure is that the other authors go on to develop a theory of discrete exterior calculus of forms only by 
introducing interpolation of forms, which we will be able to avoid. The formal definition of discrete forms 
follows. 

Definition 5.2. A primal discrete k-form a is a homomorphism from the chain group Ck{K;1) to the 
additive group R. Thus, a discrete k-form is an element of Hom(Ck(K), R), the space of cochains. This 
space becomes an abelian group if we add two homomorphisms by adding their values in R. The standard 
notation for Hoia(Ck(K),R) in algebraic topology is C k (K;R). But we will often use the notation Vt\{K) 
for this space as a reminder that this is the space of discrete (hence the d subscript) k-forms on the simplicial 
complex K . Thus, 

fl k d (K) := C k (K;R) = Rom(C k (K),R) . 
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Note that, by the above definition, given a fc-chain J^i a i c i (where a, G Z) and a discrete fc-form a, we 
have that 

and for two discrete fc-forms a, (3 G f^(if) and a fc-chain c G Ck(K; Z), 

(a + /3)(c) = a(c) + /3(c) . 

In the usual exterior calculus on smooth manifolds integration of fc-forms on a fc-dimensional manifold is 
defined in terms of the familiar integration in M. k . This is done roughly speaking by doing the integration in 
local coordinates, and showing that the value is independent of the choice of coordin ates, due to the chang e 
of variables theorem in M fc . For details on this, see the first few pages of Chapter 7 of lAbraham et all 1988]. 
We will not try to introduce the notion of integration of discrete forms on a simplicial complex. Instead the 
fundamental quantity that we will work with is the natural bilinear pairing of cochains and chains, defined 
by evaluation. More formally, we have the following definition. 

Definition 5.3. The natural pairing of a k-form a and a k-chain c is defined as the bilinear pairing 

(a, c) = a(c). 

As mentioned above, in discrete exterior calculus, this natural pairing plays the role that integration of 
forms on chains plays in the usual exterior calculus on smooth manifolds. The two are related by a procedure 
done at the time of discretization. Indeed, consider a simplicial triangulation K oi & polyhedron in R™, i.e., 
consider a "flat" discrete manifold. If we are discretizing a continuous problem, we will have some smooth 
forms defined in the space \K\ C R n . Consider such a smooth fc-form a k . In order to define the discrete 
form corresponding to a k , one would integrate a k on all the fc-simplices in K. Then, the evaluation of a 1 ^ 
on a fc-simplex a k is defined by a k t (cr k ) := J ak a k . Thus, discretization is the only place where integration 
plays a role in our discrete exterior calculus. 

In the case of a non-flat manifold, the situation is somewhat complicated by the fact that the smooth 
manifold, and the simplicial complex, as geometric sets embedded in the ambient space do not coincide. A 
smooth differential form on the manifold can be discretized into the cochain representation by identifying 
the vertices of the simplicial complex with points on the manifold, and then using a local chart to identify 
fc-simplices with A:- volumes on the manifold. 

There is the possibility of fc-volumes overlapping even when their corresponding fc-simplices do not inter- 
sect, and this introduces a discretization error that scales like the mesh size. One can alternatively construct 
geodesic boundary surfaces in an inductive fashion, which yields a partition of the manifold, but this can be 
computationally prohibitive to compute. 

Now we can define the discrete exterior derivative which we will call d, as in the usual exterior calculus. 
The discrete exterior derivative will be defined as the dual, with respect to the natural pairing defined above, 
of the boundary operator, which is defined below. 

Definition 5.4. The boundary operator dk '■ C'k(K;Z) — > Ck—i{K',Z) is a homomorphism defined by its 
action on a simplex a k = [vo, . . . , Vk], 

k 

dk<J k = d k ([v , . . .,v k }) = ^(-1)> , ...,Vi,...,Vk], 

i=0 

where [vq, . . . ,Vi, . . . , Vk] is the (fc — I) -simplex obtained by omitting the vertex Vi. Note that dk ° dk+i = 0. 

Example 5.2. Given an oriented triangle [vq, Vi, 1)2] the boundary, by the above definition, is the chain 
[vifVz] — [vo,V2] + [vo,vi], which are the three boundary edges of the triangle. 
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Definition 5.5. On a simplicial complex of dimension n, a chain complex is a collection of chain groups 
and homomorphisms d k , such that, 

> C n (K) -^-> ••• C k (K) ••• -*-> C (X) > 0, 

and <9 fe o <9 fe+1 = 0. 

Definition 5.6. The coboundary operator. S k : C k (K) — ► C is defined by duality to the boundary 

operator, with respect to the natural bilinear pairing between discrete forms and chains. Specifically, for a 
discrete form a k G O^(if), and a chain Ck+i S Ck+i{K;'L), we define 8 k by 

(5.1) (8 k a k ,c k+1 ) = (a k ,d k+1 c k+1 ) . 

That is to say 

6 k (a k ) = a k od k+1 . 
This definition of the coboundary operator induces the cochain complex, 



< C n {K) ^— ■■■ C k (K) <r-^ ■■■ C°(K) i 0, 

where it is easy to see that S k+1 o S k = 0. 

Definition 5.7. The discrete exterior derivative denoted by d : Qh(K) — > r^ +1 (i4T) is defined to be the 
coboundary operator 5 k . 

Remark 5.1. With the above definition of the exterior derivative, d : il^K) — > fl^ 1 ^), and the rela- 
tionship between the natural pairing and integration, one can regard eauation \5.1\ as a discrete generalized 
Stokes' theorem. Thus, given a k-chain c, and a discrete k-form a, the discrete Stokes' theorem, which is 
true by definition, states that 

(da, c) = {a, dc) . 
Furthermore, it also follows immediately that d fc+1 d fe = 0. 

Dual Discrete Forms. Everything we have said above in terms of simplices and the simplicial complex K 
can be said in terms of the cells that are duals of simplices and elements of the dual complex *K . One just 
has to be a little more careful in the definition of the boundary operator, and the definition we construct 
below is well-defined on the dual cell complex. This gives us the notion of cochains of cells in the dual 
complex and these are the dual discrete forms. 

Definition 5.8. The dual boundary operator, d k : C k (*K; Z) — > C k -\ (*K; Z), is a homomorphism defined 
by its action on a dual cell & k = -ka n ~ k = *[vo, ■ ■ ■ , v n - k ], 

da k =d*{v , Vn-k] 

where - n ^ k + 1 i s oriented so that it is consistent with the induced orientation on a n ~ k . 

6. Hodge Star and Codifferential 

In the exterior calculus for smooth manifolds, the Hodge star, denoted *, is an isomorphism between the 
space of fc-forms and (n — fc)-forms. The Hodge star is useful in defining the adjoint of the exterior derivative 
and this is adjoint is called the codifferential. The Hodge star, * : il k (M) — > tt n ~ k (M), is in the smooth 
case uniquely defined by the identity, 

((a k ,l3 k ))v = a k A*/3 fe , 
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where (( , )) is a metric o n differential forms, an d v is the volume-form. For a more in-depth discussion, see, 
for example, page 411 ofllb raham et 

id mm- 

The appearance of k and (n — k) in the definition of Hodge star may be taken to be a hint that primal 
and dual meshes will play some role in the definition of a discrete Hodge star, since the dual of a /c-simplex 
is an (n — k)-ce\\. Indeed, this is the case. 

Definition 6.1. The discrete Hodge Star is a map * : Q!^(K) — > tt d ~ k (-kK), defined by its action on 
simplices. For a k-simplex o~ k , and a discrete k-form a k , 

T \ 1 {a k ,a k ) = : 1 , . (*« fc ,*cr fc ). 
|cr fe p ' \-ka k y ' 

The idea that the disc rete Hodge star maps primal discrete forms to dual forms, and vice versa, is well- 
known. See, for example. ISen et all 200fll | . However, notice we now make use of the volume of these primal 



and dual meshes. But the definition we have given above does appear in the work o f lHiptmairl |2002j 

The definition implies that the primal and dual averages must be equal. This idea has already been 
introduced, n ot in the context of e xterior calculus, but in an attempt at defining discrete differential geometry 
operators, see lMever et alJ 2002j . 



Remark 6.1. Although we have defined the discrete Hodge star above, we will show in Remark \12.1\ of J7 
that if an appropriate discrete wedge product and metric on discrete k-forms is defined, then the expression 
for the discrete Hodge star operator follows from the smooth definition. 

Lemma 6.1. For a k-form a k , 

, , k (a \k(n — k) k 

* * a = {— 1) v 'a . 

Proof. The proof is a simple calculation using the property that for a simplex or a cell a k , (er fe ) = 
^_ 1 )fe(n-fe) cr fc (Equation ED . □ 

Definition 6.2. Given a simplicial or a dual cell complex K the discrete codifferential operator, S : 

tt k+1 (K) -> n k (K), is defined by 5{n° d (K)) = and on discrete (k + 1) -forms to be 

5/3 = {-l) nk+1 *d*[3. 

With the discrete forms, Hodge star, d and 6 defined so far, we already have enough to do an interesting 
calculation involving the Laplace-Beltrami operator. But, we will show this calculation in Rafter we have 
introduced discrete divergence operator. 

7. Maps between 1-Forms and Vector Fields 

Just as discrete forms come in two flavors, primal and dual (being linear functionals on primal chains 
or chains made up of dual cells), discrete vector fields also come in two flavors. Before formally defining 
primal and dual discrete vector fields, consider the examples illustrated in Figure The distinction lies in 
the choice of basepoints, be they primal or dual vertices, to which we assign vectors. 

Definition 7.1. Let K be a flat simplicial complex, that is, the dimension of K is the same as that of the 
embedding space. A primal discrete vector field X on a flat simplicial complex K is a map from the 
zero- dimensional primal subcomplex (i.e., the primal vertices) to M. N . We will denote the space of such 
vector fields by Xd(K). The value of such a vector field is piecewise constant on the dual n-cells of *K . 
Thus, we could just as well have called such vector fields dual and defined them as functions on the n-cells 
of*K. 

Definition 7.2. A dual discrete vector field X on a simplicial complex K is a map from the zero- 
dimensional dual subcomplex (*K) ^ ( i. e, the circumcenters of the primal n simplices ) to S, N such that its 
value on each dual vertex is tangential to the corresponding primal n-simplex. We will denote the space of 
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Figure 4. Discrete vector fields. 



such vector fields by Xd{*K). The value of such a vector field is piecewise constant on the n-simplices of 
K . Thus, we could just as well have called such vector fields primal and defined them as functions on the 
n-simplices of K . 

Remark 7.1. In this paper we have defined the primal vector fields only for flat meshes. We will address 
the issue of non-flat meshes in separate work. 

As in the smooth exterior calculus, we want to define the flat (b) and sharp (ft) operators that relate forms 
to vector fields. This allows one to write various vector calculus identities in terms of exterior calculus. 

Definition 7.3. Given a simplicial complex K of dimension n, the discrete flat operator on a dual 
vector field, b : Xd{*K) —> Q d (K), is defined by its evaluation on a primal 1 simplex a 1 , 

where X ■ a 1 is the usual dot product of vectors in M. N , and a 1 stands for the vector corresponding to a 1 , 
and with the same orientation. The sum is over all a n containing the edge a 1 . The volume factors are in 
dimension n. 

Definition 7.4. Given a simplicial complex K of dimension n, the discrete sharp operator on a primal 
1-form, j) : VL d (K) — > Xd{*K), is defined by its evaluation on a given vertex v, 

Jiv) = ^2 (a,[v,a }) - n [v,<T°] , 

[v,a°] (T n y[v,(T ] ' 

where the outer sum is over all l-simplices containing the vertex v, and the inner sum is over all n-simplices 
containing the 1-simplex [v,a ]. The volume factors are in dimension n, and the vector fWo-oi is the normal 
vector to the simplex [v,a a ], pointing into the n-simplex a n . 

For a discussion of the proliferation of discrete sharp and flat operators that arise from considering the 
interpolation of differential forms and vector fields, please see iHiranl |2003| . 
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8. Wedge Product 

As in the smooth case, the wedge product we will construct is a way to build higher degree f orms from lower 
degree ones. For information about the smooth case, see the first few pages of Chapter 6 of lAbraham et afl 
19881. 



Definition 8.1. Given a primal discrete k-form a € Qg(K), and a primal discrete l-form f3 l E ft l d (K), the 
discrete primal-primal wedge product, A : flH(K) x fl l JK) — > fl d +l (K), defined by the evaluation on a 
(k + I) -simplex a k+l — [vq, ■ ■ ■ , Vk+l] is given by 



T k+l 



where Sk+i+i is the permutation group, and its elements are thought of as permutations of the numbers 
0,...,k + 1 + 1. The notation r(a k+l ) stands for the simplex [u r (o) 5 ■ • ■ ,v T (k+i)}- Finally, the notation 
a ■-— ■ /3(r(cr fc+i )) is borrowed from algebraic topology (see, for example, page 206 of \Hatche~\ \200i\l ) and is 
defined as 

a — f3{T(cr k+l )) := (a, [v t(0) , . . . , v r(k) ] }(/?, [w T (k)> . . . , v r{k+ i)] ) . 

Example 8.1. When we take the wedge product of two discrete 1- forms, we obtain terms in the sum that 
are graphically represented in Figure [3| 



Figure 5. Terms in the wedge product of two discrete 1-forms. 



Definition 8.2. Given a dual discrete k-form a k G QH(*K), and a primal discrete l-form f3 l 6 Vl l d {-kK), the 
discrete dual-dual wedge product, A : Q d (*K) x tt l d (*K) — > fl d +l (-kK), defined by the evaluation on a 
(k + l)-cell a k+l = -ka n ^ k ~ l , is given by 

(a k A/3 l ,a k+l ) ={a k A (3 l ,*o n ~ k ~ l ) 

= X! si S n ( cr ™~ fc ~ i ! i v k+l, ■ ■ ■ ,v n ]) ^2 sign(r) 

■ (a fc ,*[u r ( ), . . . ,v T{ i_ 1) ,v k +u ■ ■ ■ ,Vn]}0 l ,*[v T {i), ■ ■ ■ ,v T ( k+l ^x),Vk+i, ■ ■ -,v n ]) 
where o~ n — [vq, ■ ■ ■ ,V n ], and, without loss of generality, assumed that <j n ~ k ~ l = ±[v k +i, . . . ,V n ], 

Anti-Commutativity of the Wedge Product. 

Lemma 8.1. The discrete wedge product, A : C k (K) x C l (K) — > C k+l {K), is anti- commutative, i.e., 

a k A(3 k = {-l) kl l3 l Aa k . 



18 



MATHIEU DESBRUN, ANIL N. HIRANI, MELVIN LEOK, AND JERROLD E. MARSDEN 



Proof. We first rewrite the expression for the discrete wedge product using the following computation, 
E sign(f)|cr fc+z n*^ (fe) |(a fc ,f[wo, . . . , v k ])(3 l , f [v k , . . .,v k +i\) 

= E {~ 1 ) k ~ 1 sign(f)\a k+l n*v f{k) \ {a k ,f[vi,...,vo,Vk])(J?,f[v k ,...,Vk+i\) 

= E ( _1 ) fe_1 sign(f)|cr fe+ ' n *v fp( p ) \(a k ,fp[v 1 , . . . , v k ,v ])((3 l ,fp[v , v k+1 , v k+ i}) 

= E (-l) fc_1 (-l) fc sign(f)|cr fe +' n *v fp(0) \(a k ,Tp[v , v k ])(f3 l ,fp[v , v k+1 , . . . , v k+l \) 

- E {-i) k -\-i)\-i)^{fp)W k+l n^ (0) | 

rpGSk+i + ip 

■ {a k ,fp[v , . . .,v k ])(P l ,fp[v ,v k+ i, . . . ,v k+ i}} 
= E siga(T)\a k+l n *v T(0) \(a k ,T[v , . . . , v k ])(/3 l , t[v , v k+1 , . . .,v k+l }) . 

Here, we used the elementary fact, from permutation group theory, that a k + 1 cycle can be written as 
the product of k transpositions, which accounts for the (— l) fe factors. Also, p is a transposition of and k. 
Then, the discrete wedge product can be rewritten as 



(a k A f3 l , a k+l ) = Tj~^jy^ E W^';*^ 1 

• (a k , [v r (o), ■ ■ ■ ,v T ( k) ]){[3\ [w T ( ),r(fe+i), • ■ • ,v T(k+l) ]). 
For ease of notation, we denote [v , ■ ■ ■ , v k ] by a k , and [v , v k+i , . . . , v k+ i] by a 1 . Then, we have 

<« fc A/3V+') = -l- E sign(r)^^(a fe ,r(a fe ))(/3',V)). 

Furthermore, we denote [u , Vi + i . . . , v k+ {\ by a k , and [«o, v\, . . . , vi] by a 1 . Then, 

(P l A a k , a k+l ) = 77~]\~\ E Bign(f) |gfc+; ^ (0)l (a fc ,f(a fc ))(/3',f(a')). 

Consider the permutation G 5fe + ; + i, given by 

'0 1 ... k k + 1 ... k + l^ 
v l + l ... k + l 1 

which has the property that 

a k = 0(a k ), 
a 1 = 6{a l ). 

Then, we have 

{/3 l A a k , a k+l ) = Tj~^~jT\ E «W '^T^ ' <"*■ rfrW.ffr 1 )) 

= TFTTV E si g n(.) ^^ ^,^(^))^^(.')) 
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£ sign(f0) sign(fl) |gfc+ ' n ^ (0)l (a k , ?9{a k )> ( P l ,?9(a% 



(k + l)) 

y ' fees k+l+1 e 

By making the substitution, r = f9, and noting that Sk+i+iO = Sk+i+i, we obtain 

(^Aa fc ,^>= S ign(^)-l- J2 sign(r)^-^|^(a fe ,r( ( T fc ))(/3',r(a i )} 

TGSfc + 1 + 1 

= sign(#)(a fc A f3 l , a h+l ) . 
To obtain the desired result, we simply need to compute the sign of 6, which is given by 

sign(0) = (-1)*'. 

This follows from the observation that in order to move each of the last I vertices of a k+l forward, we require 
k transpositions with V\, . . . ,Vk- Therefore, we obtain 

(f3 l Aa k ,a k+l ) =sign(0)(a fc A/3 l ,<j k+l ) = (-l) kl {a k A (3\ a k+l ), 

and 

a k Af3 l = {-l) kl f3 l Aa k . □ 



Leibniz Rule for the Wedge Product. 

Lemma 8.2. The discrete wedge product satisfies the Leibniz rule, 

d(a k A p l ) = (da k ) A ft + (-l) k a k A (d(3 l ). 

Proof. The proof of the Leibniz rule for discrete wedge products is directly analogous to the pro of of the 
cobou ndary formula for the simplicial cup product on cochains, which can be found on page 206 of lHatcherl 
2001J. This is because the discrete exterior derivative is precisely the coboundary operator, and the wedge 
product is constructed out of weighted sums of cup products. 

The cup product satisfies the Leibniz rule for an given partial ordering of the vertices, and the permutations 
in the signed sum in the discrete wedge product correspond to different choices of partial ordering. We then 
obtain the Leibniz rule for the discrete wedge product by applying it term- wise for each choice of permutation. 

Consider 

fe+i 



i=0 K >' reS t+1 + i 1 1 



and 



{a k , [u t (q), ...,Vi,... , i> T (fc+i)] )(/?', [iVffc+i), . . . ,« T (/s + ; + i)]), 



(-i) k {a k A(df3 l ),„ k+l+1 ) = (-i) k £ i-^jr- ni E si § n ( r ) 



(k + l)\ ^ \o- k+l \ 

■ (a k , [v T (o), ■ ■ ■,v T (k)})(l3 l , [«r(fc)j ...,Vi,.. . ,v T (fc+i+i)])- 
The last set of terms, i = k + 1, of the first expression cancels the first set of terms, i = k, of the second 
expression, and what remains is simply (a k A (3 l , do~ k+l+1 ). Therefore, we can conclude that 

((da fc ) A [3\ a k+l+1 ) + (-l) k (a k A {dl3 l ),a k+l+1 ) = (a k A /3 l , d<j k+l+1 ) = (d(a k A (3 l ), <j k+l+1 ), 
or simply that the Leibniz rule for discrete differential forms holds, 

d(a fe A f3 l ) = (da k ) A j3 l + {-l) k a k A (d/3 z ). □ 
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Associativity for the Wedge Product. The discrete wedge product which we have introduced is not 
associative in general. This is a consequence of the fact that the stencil for the two possible triple wedge 
products are not the same. In the expression for (a k A(f3 l Aj m ),a k+l+m ), each term in the double summation 
consists of a geometric factor multiplied by (a , <j k )((3 l , a l )(j m , <j m ) for some k, I, m simplices er fe , a 1 , a m . 

Since [3 l and 7™ 1 are wedged together first, a 1 and a m will always share a common vertex, but a k could 
have a vertex in common with only a , or only <j m , or both. We can represent this in a graph, where the 
nodes denote the three simplices, which are connected by an edge if, and only if, they share a common vertex. 
The graphical representation of the terms which arise in the two possible triple wedge products are given in 
Figure 




a A (/? A 7) {a A (f) A 7 



Figure 6. Stencils arising in the double summation for the two triple wedge products. 

For the wedge product to be associative for all forms, the two stencils must agree. Since the stencils for 
the two possible triple wedge products differ, the wedge product is not associative in general. However, in 
the case of closed forms, we can rewrite the terms in the sum so that all the discrete forms are evaluated on 
triples of simplices that share a common vertex. This is illustrated graphically in Figure [7] 




FIGURE 7. Associativity for closed forms. 

This result is proved rigorously in the follow lemma. 

Lemma 8.3. The discrete wedge product is associative for closed forms. That is to say, for a k £ C k (K), 
P l £ C l {K), 7™ £ C m (K), such that da k = 0, d[3 l = 0, d 7 m = 0, we have that 

(a k A/3') A7" 1 = a k A (f3 l A7" 1 ). 

Proof. 

((a k A/3 ; ) A7 m ,(j fc+ ' +m ) 
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Y sign(T)(a k A(3 l ,T[v ,...,v k+ i])(^ m ,T[v k+ i,...,v k+ i +m }) 

= Y Y sign(T)sign(p)(a fc ,,OT[t>o, . . . ,v k ]} 

TeSfc+i+m+i peSfe+i+i 

• (P l ,pr[vk, ■ ■ ■ ,v k+ i]}(j m ,T[v k+ i, . . .,v k+l+m ]) 

Here, either pr(k) = r(k + I), in which case all three permuted simplices share v T ( k +i) as a common vertex, 
or we need to rewrite either {a k ,pr[v , ■ ■ ■ ,v k ]) or {[3 l ,pr[v k , . . . ,v k+ i]), using the fact that a k and (3 l are 
closed forms. 

If v r(k+i) ^ pr[vo, ■ ■ ■ ,v k ], then we need to rewrite (a k , pr[vo, . . . , v k ]) by considering the simplex obtained 
by adding the vertex v T ^+i) to pr[vo, . . . , v k ], which is [v T (k+i),v pT (o), . . . , v pT (k)]- Then, since a k is closed, 
we have that 

= (da k , K( fe+ /),Wp T (o), • • • ,v P r(k)}) 
= (a fe ,9K(fe+o> w pr(o), • • • ,v pT (k)}) 

k 

= (a k , bpr(O), • ■ • ,V pr (k)}) - J^C - !)'^! [ V T(k+l),V pT (0), ■ ■ ■ ,V pr (i), ■ ■ ■ ,Vpr(k)]) 

or equivalently, 

k 

(« fe , bpr(O), • ■ •, Vpr(k)}) = X]^ -1 )'^' K(fe+0> W pr(0), • ■ •, V(i)' ' ' -' V pr(k)])- 

i=0 

Notice that all the simplices in the sum, with the exception of the last one, will share two vertices, v T (fc+i) 
and ^ pr (A;) with pr[v k , . . . , v k+l ], and so their contribution in the triple wedge product will vanish due to the 
anti-symmetrized sum. 

Similarly, if v T (fc+i) ^ P T i v k, ■ ■ ■ ,Vk+l]> using the fact that /3 is closed yields 

k+i 

(a k , [Vpr(k), ■ ■ ■ , Vpr(k+l)}) = Y^ 1 ^^ ^ > i V r(k+l),Vpr(k),- ■ ■, V(i)' ' ' ' ' V pr(k+l)})- 

i=k 

As before, all the simplices in the sum, with the exception of the last one, will share two vertices, v T (fc+i) 
and Wp T (fe) with pr[v , ■ ■ . , v k ], and so their contribution in the triple wedge product will vanish due to the 
anti-symmetrized sum. 

This allows us to rewrite the triple wedge product in the case of closed forms as 

k-\-l+m 

((« t A/3 i )A7 m ,^ +I+ '")= Y Y sign(p i r)(a k ,p i r[v ,...,Vk})(p l ,p i T[v ,v k+ i,...,Vk + i}) 

i=0 reS k+ i + m 

■ il" 1 , PiT[vQ,V k +l+l, ■ ■ ■ ,Vk+l+m]) , 

where r G S k+ i +m is thought of as acting on the set {1, . . . , k + I + m}, and pi is a transposition of and i. 
A similar argument allows us to write a k A ((3 l A 7" 1 ) in the same form, and therefore, the wedge product is 
associative for closed forms. □ 



Remark 8.1. This lemma is significant, since if we think of a constant smooth differential form, and 
discretize it to obtain a discrete differential form, this discrete form will be closed. As such, this lemma 
states that in the infinitesimal limit, the discrete wedge product we have defined will be associative. 
In practice, if we have a mesh with characteristic length Ax, then we will have that 

1 k t l+m M k A (f3 l A 7 m ) - (a k A/3') A 7 "> fc + (+m > - O(Ax), 
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which is to say that the average of the associativity defect is of the order of the mesh size, and therefore 
vanishes in the infinitesimal limit. 

9. Divergence and Laplace-Beltrami 

In this section, we will illustrate the application of some of the DEC operations we have previously defined 
to the construction of new discrete operators such as the divergence and Laplace-Beltrami operators. 

Divergence. The divergence of a vector field is given in terms of the Lie derivative of the volume-form, 
by the expression, (div(X)/i = £xM- Physically, this corresponds to the net flow per unit volume of an 
infinitesimal volume about a point. 

We will define the discrete divergence by using the formulas defining them in the smooth exterior calculus. 
The divergence definition will be valid for arbitrary dimensions. The resulting expressions involve operators 
that we have already defined and so we can actually perform some calculations to express these quantities 
in terms of geometric quantities. We will show that the resulting expr ession in terms of geometric quantities 
is the same as that derived by variational means in lTong et alJ|2003| . 

Definition 9.1. For a discrete dual vector field X the divergence div(A) is defined to be 

div(X) = -SX h . 

Remark 9.1. The a bove definition is a theorem in smooth exterior calculus. See, for example, page 458 of 
\Abraham et~ai\ \l98&l . 

As an example, we will now compute the divergence of a discrete dual vector field on a two-dimensional 
simplicial complex K, as illustrated in Figure |H1 




Figure 8. Divergence of a discrete dual vector field. 



A similar derivation works in higher dimensions, where one needs to be mindful of the sign that arises 
from applying the Hodge star twice, * *a k = {— l) fc (™~ fc )a fe . Since div(X) = — SX b , it follows that div(X) 
• (1 * X . Since this is a primal 0-form it can be evaluated on a 0-simplex er°, and we have that 

(div(a;),a } = <*d * X\ a ) . 

Using the definition of discrete Hodge star, and the discrete generalized Stokes' theorem, we get 

^<div(X),a°) = ^(**dH<A b ,^) 
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(*X\d{*a )} 



The second equality is obtained by applying the definition of the Hodge star, and the last equality is obtained 
by applying the discrete generalized Stokes' theorem. But, 



<9(*a c 



E 



as given by the expression for the boundary of a dual cell in Eauation l5.8l Thus, 



1 \- \*^\ /v \> u 

1 n lWnO 1 



★ a 1 - 



E 



E 



★ o- 1 n cr 2 



■A- (J 



■ E E 



l^a 1 n<r 2 | 
k 1 ! 



1 o-i^cr 1 1 

This expression has the nice property that the divergence theorem holds on any dual n-chain, which, 
as a set, is a simply connected subset of \K\. Furthermore, the coefficients we computed for the discrete 
divergence operator are the unique ones for which a discrete divergence theorem holds. 

Laplace— Beltrami. The Laplace-Beltrami operator is the generalization of the Laplacian to curved spaces. 
In the smooth case the Laplac e-Beltrami operator on smooth functions is defined to be V 2 = div o curl = Sd. 
See, for example, page 459 of [Ah raha m et al . 1988J. Thus, in the smooth case, the Laplace-Beltrami on 
functions is a special case of the more general Laplace-deRham operator, A : fl k (M) — ► fl k (M), defined by 
A = dS + Sd. 

As an example, we compute A/ on a primal vertex a , where / £ fl^(K), and if is a (not necessarily 
flat) triangle mesh in R 3 , as illustrated in Figure 
This calculation is done below. 

^(Af,a°) = (Sdf.a ) 

= -(*d*df,a°) 
1 



1 



^<d*d/,*a°) 



mi/, 



~k G 



E w.- 1 ) 
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Figure 9. Laplace-Beltrami of a discrete function. 



1 1 o-Vct 1 ■ 

where 9a 1 = u — er°. But, the above is the same as the formula involving cotangents found bv lMever et alJ 
|2002j without using discrete exterior calculus. 

Another interesting aspect, which will be discussed in ^121 is that the characterization of harmonic func- 
tions as those functions which vanish when the Laplace-Beltrami operator is applied is equivalent to that 
obtained from a discrete variational principle using DEC as the means of discretizing the Lagrangian. 

10. Contraction and Lie Derivative 

In this section we will discuss some more operators that involve vector fields, namely contraction, and Lie 
derivatives. 

For contraction, we will first define the usual smooth contraction algebraically, by relating it to Hodge 
star and wedge products. This yields one potential approach to defining discrete contraction. However, since 
in the discrete theory we are only concerned with integrals of forms, we can use the interesting notion of 
extrusion of a manifold by the flow of a vector field to define the integral of a contracted discrete differential 
form. 

We learned about this definition of contraction via extrusion from lBossavitl |2002bj , who goes on to define 
discrete extrusion in his paper. Thus, he is able to obtain a definition of discrete contraction. Extrusion 
turns out to be a very nice way to define integrals of operators involving vector fields, and we will show how 
to define integrals of Lie derivatives via extrusion, which will yield discrete Lie derivatives. 

Definition 10.1. Given a manifold M, and S, a k-dimensional submanifold of M, and a vector field 
X € X(M), we call the manifold obtained by sweeping S along the flow of X for time t as the extrusion of 
S by X for time t, and denote it by E X (S). The manifold S carried by the flow for time t will be denoted 
tp x (S). 

Example 10.1. Fiaure UHl illustrates the 2-simplex that arises from the extrusion of a 1-simplex by a discrete 
vector field that is interpolated using a linear shape function. 
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Figure 10. Extrusion of 1-simplex by a discrete vector field. 



Contraction (Extrusion). We first establish an integral property of the contraction operator. 
Lemma 10.1. 



d_ 

lit 



ft 



t=0-JE' x (S) 



Proof. Prove instead that 



dr 



0. 



Then, by first fundamental theorem of calculus, the desired result will follow. To prove the above, simply 
take coordinates on S and carry th em along with th e flow and define the transversal coordinate to be the 
flow of X. This proof is sketched in lBossavfj |2002bj |. □ 

This lemma allows us to interpret contraction as being the dual, under the integration pairing between 
fc-forms and A-volumes, to the geometric operation of extrusion. The discrete contraction operator is then 
given by 

d 



<i x a 



t=0 



where the evaluation of the RHS will typically require that the discrete differential form and the discrete 
vector field are appropriately interpolated. 

Remark 10.1. Since the dynamic definition of the contraction operator only depends on the derivative of 
pairing of the differential form with the extruded region, it will only depend on the vector field in the region 
S, and not on its extension into the rest of the domain. 

In addition, if the interpolation for the discrete vector field satisfies a superposition principle, then the 
discrete contraction operator will satisfy a corresponding superposition principle. 



Contraction (Algebraic). Contraction is an operator that allows one to combine vector fields and forms. 
For a smooth manifold M, the contraction of a vector field X 6 %(M) with a (k + l)-form a £ ft k+1 (M) 
is written as ixct, and for vector fields X\, . . . ,Xk £ X(M), the contraction in smooth exterior calculus is 
defined by 

i x a(X 1 ,...,X k )=a(X,X 1 ,...,X k ). 
We define c ontraction by using an identity that is true in smooth exterior calculus. This identity originally 
appeared in iHiranil [2003], and we state it here with proof. 
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Lemma 10.2 (Hirani 2003]). Given a smooth manifold M of dimension n, a vector field X 6 3L(M), and 
a k-form a 6 fl k (M), we have that 

i z a=(-l) tM, *(*aM l ). 

Proof. Recall that for a smooth function / 6 fl°(M), we have that ixa = fix®. This, and the multilinearity 
of a, implies that it is enough to show the result in terms of basis elements. In particular, let t G S n be a 
permutation of the numbers 1, . . . n, such that t(1) < . . . < r(fc), and r(fc + 1) < . . . < r(n). Let X — e T u\, 
for some j £ 1, . . . , n. Then, we have to show that 

ie T(3 .,e T(1) A ... A e rW = (-l) fc ("- fe ) * (*( e r « A ... A e r(fc) ) A e r(j) ) . 

It is easy to see that the LHS is if j > k, and it is 

(-ir 1 (e r(1) A...A^...Ae jW ), 

otherwise, where e T ^) means that e T V' is omitted from the wedge product. Now, on the RHS of Equa- 
tion ^3 we have that 

*(e r(1) A ... A e T «) = sign(r)(e T (' £ + 1 ) A ... A e T ^) . 

Thus, the RHS is equal to 

(_!)*(»•-*) sign ( r ) * ( e r(fc+l) A _ _ A e r(n) A e r(j)) ^ 

which is as required if j > /c. So, assume that 1 < j < k. We need to compute 

*(e T(fc+1) A... Ae T <") Ae T <^), 

which is given by 

se r WA...e^)A...Ae T W, 
where the sign s = ±1, such that the equation, 

s e r( - k+1 1 A ... A e T (") A e r(j) A e T ^ A ... A e^) A ... A e T ( fe ) = /i , 
holds for the standard volume-form, /i = e 1 A . . . A e". This implies that 

fl = (-l) J '- 1 (-l)*( n - fc) sign(r). 
Then, RHS = LHS as required. □ 

Since we have expressions for the discrete Hodge star (*), wedge product (A), and flat (b), we have the 
necessary ingredients to use the algebraic expression proved in the above lemma to construct a discrete 
contraction operator. 

One has to note, however, that the wedge product is only associative for closed forms, and as a consequence, 
the Leibniz rule for the resulting contraction operator will only hold for closed forms as well. This is, however, 
sufficient to establish that the Leibniz rule for the discrete contraction will hold in the limit as the mesh is 
refined. 

Lie Derivative (Extrusion). As was the case with contraction, we will establish a integral identity that 
allows the Lie derivative to be interpreted as the dual of a geometric operation on a volume. This involves 
the flow of a volume by a vector field, and it is illustrated in the following example. 

Example 10.2. Figure 1771 illustrates the flow of a l-simplex by a discrete vector field interpolated using a 
linear shape function. 

Lemma 10.3. 
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Figure 11. Flow of a 1-simplcx by a discrete vector field. 



Proof. 



Is Jo 

ct , 



F t *(£xP) = 


dt * H 


F* T {£ x l3)dT = 


F*P - 


F* t {£ x I3)oIt = 


f F*P — 




Js 


£x/3dr = 






J V X (S) 





t 

' 



0. □ 



This lemma allows us to define a discrete Lie derivative as follows, 

d 



where, as before, evaluating the RHS will require the discrete differential form and discrete vector field to 
be appropriately interpolated. 

Lie Derivative (Algebraic). Alternatively, as we have expressions for the discrete contraction operator 
(ix), and exterior derivative (d), we can construct a discrete Lie derivative using the Cartan magic formula, 

£x^> = ixduj + dixw. 

As is the case with the algebraic definition of the discrete contraction, the discrete Lie derivative will only 
satisfy a Leibniz rule for closed forms. As before, this is sufficient to establish that the Leibniz rule will hold 
in the limit as the mesh is refined. 



11. Discrete Poincare Lemma 

In this section, we will prove the discrete Poincare lemma by const ructing a homotopy o perator though a 
generalized cocone construction. This section is based on the work in lDesbrun et al.l |2003j| . 

The standard cocone construction fails at the discrete level, since the cone of a simplex is not, in general, 
expressible as a chain in the simplicial complex. As such, the standard cocone does not necessarily map 
fc-cochains to (k — l)-cochains. 

An example of how the standard cone construction fails to map chains to chains is illustrated in Figure lT^l 
Given the simplicial complex on the left, consisting of triangles, edges and nodes, we wish, in the center 
figure, to consider the cone of the bold edge with respect to the top most node. Clearly, the resulting cone in 
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the right figure, which is shaded grey, cannot be expressed as a combination of the triangles in the original 
complex. 

AMA 

Figure 12. The cone of a simplex is, in general, not expressible as a chain. 

In this subsection, a generalized cone operator that is valid for chains is developed which has the essential 
homotopy properties to yield a discrete analogue of the Poincare lemma. 

We will first consider the case of trivially star-shaped complexes, followed by logically star-shaped com- 
plexes, before generalizing the result to contractible complexes. 

Definition 11.1. Given a k-simplex a k = [vq, ... ,1^] we construct the cone with vertex w and base a k , as 
follows, 

woa k = [w,v , . . .,v k ]. 
Lemma 11.1. The geometric cone operator satisfies the following property, 

d{woa k ) +wo(da k ) = a k . 
Proof. This is a standard result from simplicial algebraic topology. □ 

Trivially Star-Shaped Complexes. 

Definition 11.2. A complex K is called trivially star-shaped if there exists a vertex w G such that 

for all a k G K , the cone with vertex w and base o~ k is expressible as a chain in K . That is to say, 

3w E X (0) | Vcr k e K , w o a k E C k+1 (K). 

We can then denote the cone operation with respect to w as p : Ck{K) — > Ck+i(K). 

Lemma 11.2. In trivially star-shaped complexes, the cone operator, p : Ck(K) — > Ck+i(K), satisfies the 
following identity, 

pd + dp = I, 

at the level of chains. 

Proof. Follows immediately from the identity for cones, and noting that the cone is well-defined at the level 
of chains on trivially star-shaped complexes. □ 

Definition 11.3. The cocone operator, H : C k (K) — > C k ~ 1 (K), is defined by 

(Ha k ,a k - 1 ) = (a k ,p(a k - 1 )). 

This operator is well-defined on trivially star-shaped simplicial complexes. 

Lemma 11.3. The cocone operator, H : C k (K) — > C k ~ 1 (K), satisfies the following identity, 

Hd + dH = I, 

at the level of cochains. 




DISCRETE EXTERIOR CALCULUS 



29 



Proof. A simple duality argument applied to the cone identity, 

pd + dp = I, 

yields the following, 

(a k ) a k ) = (a k ,{pd + dp)a k ) 

= (a k ,pda k ) + (a k ,dpa k ) 
= (Ha k ,da k ) + (da k ,pa k ) 
= ((dHa k ,a k ) + (Hda k ,a k ) 
= ({dH + Hd)a k ,o- k ). 

Hd + dH = I, 



Therefore, 
at the level of cochains. 



□ 



Corollary 11.4 (Discrete Poincare Lemma for Trivially Star-shaped Complexes). Given a closed cochain 
a k , that is to say, da k — 0, there exists a cochain (3 k ~ 1 , such that, df3 k ~ 1 — a k . 



Proof. Applying the identity for cochains, 
we have, 



Hd + dH 



(a k ,a k ) = ((Hd + dH)a k ,a k ) , 

but, da k = 0, so, 

(a k ,a k ) = (d(Ha k ),a k ). 
Therefore, f3 k ~ 1 — Ha k is such that d/?' c ~ 1 = a k at the level of cochains. 



□ 



Example 11.1. We demonstrate the construction of the tetrahedralization of the cone of a (n — l)-simplex 
over the origin. 

If we denote by v k , the projection of the Vi vertex to the k-th concentric sphere, where the 0-th concentric 
sphere is simply the central point, then we fill up the cone [c,Vi, ...v n ] with simplices as follows, 

[vi,v\, . . [vl,v\, . . .,«*], [vl,vl,v\, . . .,«*], . . . , [vf, . . 

Since S n ~ 1 is orientable, we can use a consistent triangulation of S n ~ l and these n-cones to consistently 
triangulate B n such that the resulting triangulation is star-shaped. 

This fills up the region to the 1st concentric sphere, and we repeat the process by leapfrogging at the last 
vertex to add [vf, and continuing the construction, to fill up the annulus between the 1st and 2nd 

concentric sphere. Thus, we can keep adding concentric shells to create an arbitrarily dense triangulation of 
a n-ball about the origin. 

In three dimensions, these simplices are given by 





[c,vl,vl,vl], 



J D u l> ^2> ^3 J j 
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Putting them together, we obtain Fiaure UM 





t 7 L 











Figure 13. Triangulation of a three-dimensional cone. 



This example is significant, since we have demonstrated that for any n-dimensional ball about a point, we 
can construct a trivially star-shaped triangulation of the ball, with arbitrarily high resolution. This allows 
us to recover the smooth Poincare lemma in the limit of an infinitely fine mesh, using the discrete Poincare 
lemma for trivially star-shaped complexes. 



Logically Star-Shaped Complexes. 

Definition 11.4. A simplicial complex L is logically star-shaped if it is isomorphic, at the level of an 
abstract simplicial complex, to a trivially star-shaped complex K . 

Example 11.2. We see two simplicial complexes, in Figure [7^} which are clearly isomorphic as abstract 
simplicial complexes. 





Figure 14. Trivially star-shaped complex (left); Logically star-shaped complex (right). 



Definition 11.5. The logical cone operator p : C k (L) — > C k+1 (L) is defined by making the following 
diagram commute, 

C k {K)^^C k+1 {K) 



C k {L)—^C k+1 {L) 
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Which is to say that, given the isomorphism ip : K — ► L, we define 

p L = <P°Pk ° y~ X - 

Example 11.3. We show an example of the construction of the logical cone operator. 




f 



This definition of the logical cone operator results in identities for the cone and cocone operator that 
follow from the trivially star-shaped case, and we record the results as follows. 

Lemma 11.5. In logically star-shaped complexes, the logical cone operator satisfies the following identity, 

pd + dp = I, 

at the level of chains. 

Proof. Follows immediately by pushing forward the result for trivially star-shaped complexes using the 
isomorphism. □ 

Lemma 11.6. In logically star-shaped complexes, the logical cocone operator satisfies the following identity, 

Hd + dH = I, 

at the level of cochains. 

Proof. Follows immediately by pushing forward the result for trivially star-shaped complexes using the 
isomorphism. □ 

Similarly, we have a Discrete Poincare Lemma for logically star-shaped complexes. 



Corollary 11.7 (Discrete Poincare Lemma for Logically Star-shaped Complexes). Given a closed cochain 
a k , that is to say, da k = 0, there exists a cochain (3 1 , such that, d(3 k 



Proof. Follows from the above lemma using the proof for the trivially star-shaped case. 



□ 



Contractible Complexes. For arbitrary contractible complexes, we construct a generalized cone operator 
such that it satisfies the identity, 

pd + dp = I, 

which is the crucial property of the cone operator, from the point of view of proving the discrete Poincare 
lemma. 
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The trivial cone construction gives a clue as to how to proceed in the construction of a generalized cone 
operator. Notice that if a a k+1 is a term in p(a k ), then p(a k+1 ) = 0. This suggests how we can use the 
cone identity to inductively construct the generalized cone operator. 

To define p(a k ), we consider a k+1 >- <r k , such that, <r k+1 and a k are consistently oriented. We apply 
pd + dp to a k+1 . Then, we have 

a k+1 = p(a k ) + p(da k+1 - a k ) + dp(a k+1 ). 

If we set p(a k+1 ) = 0, 

a k+1 = p(a k ) + p(da k+1 - a k ) + 0(0) 
= p(a k ) + p(Oa k+1 -a k ). 

Rearranging, we have 

p(a k )=a k+1 -p(da k+1 -a k ), 

and 

P (<7 k+1 ) = 0. 

We are done, so long as the simplices in the chain da k+1 — a k already have p defined on it. This then reduces 
to enumerating the simplices in such a way that in the right hand side of the equation, we never evoke terms 
that are undefined. 

We now introduce a method of augmenting a complex so that the enumeration condition is always satisfied. 

Definition 11.6. Given a n-complex K , consider a (n — l)-chain c n -i that is contained on the boundary of 
K, and is included in the one-ring of some vertex on OK. Then, the one-ring cone augmentation of K 

is the complex obtained by adding the n-cone u>oc n _i, and all its faces to the complex. 

Definition 11.7. A complex is generalized star-shaped if it can be constructed by repeatedly applying the 
one-ring augmentation procedure. 

We will explicitly show in Examples 111.41 and 111.71 how to enumerate the vertices in two and three 
dimensions. And in Examples 1 11. 61 and 111.81 we will introduce regular triangulations of R 2 and M 3 that can 
be constructed by inductive one-ring cone augmentation. 

Remark 11.1. Notice that a non-contractible complex cannot be constructed by inductive one-ring cone 
augmentation, since it will involve adding a cone to a vertex that has two disjoint base chains. This prevents 
us from enumerating the simplices in such as way that all the terms in do- k+1 — o~ k have had p defined on 
them, and we see in Example Ml. {A how this causes the cone identity, and hence the discrete Poincare lemma 
to break. 

Example 11.4. In two dimensions, the one-ring condition implies that the base of the cone consists of either 
one or two 1-simplices. To aid in visualization, consider Figure [731 




Figure 15. One-ring cone augmentation of a complex in two dimensions. 
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In the case of one \- simplex, [vq,v\\, when we augment using the cone construction with the new vertex 
w, we define, 

p([w\) = [v ,w] +p([v ]), p([v ,w]) = 0, 

p([vi,w]) = [v 0) vi,w] -p([v 0) vi]), p([v ,vi,w]) = 0. 

In the case of two 1-simplices, [vq,vi], [vq,V2], we have, 

KM) = [vo,w] +p(bo]), p([vq,w]) = 0, 

p([vi,w]) = [v ,vi,w] -p([v ,vi]), p([v ,vi, w}) = 0, 

P([v2,w]) = [v ,v 2 ,w] - p([v 0) v 2 ]), p([v ,v 2 ,w}) = 0. 

Example 11.5. We will now explicitly utilize the one-ring cone augmentation procedure to compute the 
generalized cone operator for part of a regular two-dimensional triangulation that is not logically star-shaped. 

As a preliminary, we shall consider a logically star-shaped complex, and augment with a new vertex, as 
seen in Figure TTh\ 



Figure 16. Logically star-shaped complex augmented by cone. 



We use the logical cone operator for the subcomplex that is logically star-shaped, and the augmentation 
rules in the example above for the newly introduced simplices. This yields, 
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Example 11.6. Clearly, the regular two-dimensional triangulation can be obtained by the successive appli- 
cation of the one-ring cone augmentation procedure, as the following sequence illustrates, 




which means that the discrete Poincare lemma can be extended to the entire regular triangulation of the 
plane. 



Example 11.7. We consider the case of augmentation in three dimensions. Denote by vq, the center of 
the one-ring on the two-surface, to which we are augmenting the new vertex w. The other vertices of the 
one-ring are enumerated in order, v\, . . . , v m . To aid in visualization, consider Figure \T1[ 

If the one-ring does not go completely around vq, we shall denote the missing term by [vq,vi, v m ]. The 
generalized cone operators are given as follows. 

k=0, 



KM) = i v o,w] +p([vo]), 



p([v ,w)) = 0, 
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k=L 



k=2, 



Figure 17. One-ring cone augmentation of a complex in three dimensions. 



P([vi,w}) = [vo>vi,w] -p([v ,vt]), p([v 0) vi,w]) = 0, 

p([v m ,w]) = [v ,v m ,w] -p([v ,v m ]), p([v ,v m ,w]) = 0, 



p([v 1 ,v 2 ,w}) = [v ,v 1 ,v 2 ,w] +p([v ,v 1 ,v 2 ]), p([vi,v 2 ,w]) = 0, 

p([v m -i,v m , w}) = [v ,v m -i,v m ,w] +p([v ,v m -i,v m }), p([v m ^i,v m ,w]) = 0. 

If it does go around completely, 

p([v m ,vi,w]) = [v ,v m ,vi,w] +p([v ,v m , vi]), p([v ,v m ,v-L,w]) = 0. 

Example 11.8. We provide a tetrahedralization of the unit cube that can be tiled to yield a regular tetrahe- 
dralization o/R 3 . The 3-simplices are as follows, 

[vooo ,vqoi, vow ,vw], Noi , vow , «ioo , v ioi] , Noi , vqw , v on , v W i] , 

[vow , v wo , "101 , viw] , [vow, Vqh, V W i , Vno] , [vqh , Uioi , Duo, um]. 

The tetrahedralization of the unit cube can be seen in Figure \17A 

Since this regular tetrahedralization can be constructed by the successive application of the one-ring cone 
augmentation procedure, the Discrete Poincare lemma can be extended to the entire regular tetrahedralization 
o/R 3 . 

In higher dimensions, we can extend the construction of the generalized cone operator inductively using 
the one-ring cone augmentation by choosing an appropriate enumeration of the base chain. Topologically, 
the base chain will be the cone of S n ~ 2 (with possibly an open (n — 2)-ball removed) with respect to the 
central point. 

By spiraling around S n ~ 2 , starting from around the boundary of the n — 2 ball, and covering the rest of 
S n ~ 2 , as in Figure IT§1 we obtain the higher-dimensional generalization of the procedure we have taken in 
Examples HOI and lTT"7l 

Notice that n — 2 is distinguished, since S 2 ~ 2 = S° is disjoint, which is why in the two-dimensional case, 
we were not able to use the spiraling technique to enumerate the simplices. 

Since we have constructed the generalized cone operator such that the cone identity holds, we have, 

Lemma 11.8. In generalized star-shaped complexes, the generalized cone operator satisfies the following 
identity, 

pd + dp = I, 
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at the level of chains. 

Proof. By construction. □ 

Lemma 11.9. In generalized star-shaped complexes, the generalized cocone operator satisfies the following 
identity, 

Hd + dH = I, 

at the level of cochains. 

Proof. Follows immediately from applying the proof in the trivially star-shaped case, and using the identity 
in the previous lemma. □ 

Similarly, we have a discrete Poincare lemma for generalized star-shaped complexes. 

Corollary 11.10 (Discrete Poincare Lemma for Generalized Star-shaped Complexes). Given a closed 
cochain a k , that is to say, da k — 0, there exists a cochain (3 k ~ 1 , such that, d^^ 1 = a k . 

Proof. Follows from the above lemma using the proof for the trivially star-shaped case. □ 

Example 11.9. We will consider an example of how the Poincare lemma fails in the case when the complex 
is not contractible. Consider the following trivially star-shaped complex, and augment by one vertex so as to 
make the region non- contractible, as show in Figure XWk 



(a) Trivially star-shaped complex (b) Non-contractible complex 

Figure 20. Counter-example for the discrete Poincare lemma for a non-contractible complex. 
Now we attempt to verify the identity, 

pd + dp = I, 

and we will see how this is only true up to a chain that is homotopic to the inner boundary. 
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Since the second term cannot be expressed as the boundary of a 2- chain, it will contribute a non-trivial effect, 
even on closed discrete forms, and therefore the discrete Poincare lemma does not hold for non-contractible 
complexes, as expected. 



12. Discrete Variational Mechanics and DEC 

We recall that discrete variational mechanics is based on a discrete analogue of Hamilton's principle, and 
they yield the discrete Euler-Lagrange equations. A particularly interesting property of DEC arises when it 
is used to construct the discrete Lagrangian for harmonic functions, and Maxwell's equations. 

In particular, for these examples, the following diagram commutes, 

Lagrangian dec Discrete Lagrangian 

L:TQ -> R > L d :QxQ^R 



Euler Lagrange dec Discrete Euler Lagrange 

££ : T 2 Q -> T*Q > ££ d : Q 3 -> T*Q 

Which is to say that directly discretizing the differential equations for harmonic functions, and Maxwell's 
equations using DEC results in the same expressions as the discrete Euler-Lagrange equations associated 
with a discrete Lagrangian which is discretized from the corresponding continuous Lagrangian by using DEC 
as the discretization scheme. 

This is significant, since it implies that when DEC is used to discretize these equations, the corresponding 
numerical scheme which is obtained is variational, and consequently exhibits excellent structure-preserving 
properties. 

In the variational principles for both harmonic functions and Maxwell's equations, we require the L 2 norm 
obtained from the L 2 inner product on fl k (M), which is given by 

(a k ,(3 k ) = [ aA*[3. 

The discrete analogue of this requires a primal-dual wedge product, which is given below for forms of 
complementary dimension. 

Definition 12.1. Given a primal discrete k-form a k S il^K), and a dual discrete (n — k)-form (3 n ~ k G 
Vt d l ~ k {-kK), the discrete primal-dual wedge product is defined as follows, 

(a k A$ n ~ k ,V ak ) = ^ {a k ,* k )0 n - k ,*a k ) 

\0~ K \ \ •kO~ K \ 
1 



k „k \ i Qti—k _, 

k 



-(aV fc )(/?"-W), 



where V a k is the n-dimensional support volume obtained by taking the convex hull of the simplex a and its 
dual cell -ko k . 

The corresponding L 2 inner product is as follows. 
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Definition 12.2. Given two primal discrete k-forms, a k ,(3 k £ Qh{K), their discrete L 2 inner product, 

(a k ,(3 k )d, is given by 

(^'^ = E | J^L (a k ,a k ){*[3,*o k ) 

= - T (a k ,<j k )(*f3,*o- k ). 
n 

a k eK 



Remark 12.1. Notice that it would have been quite natural from the smooth theory to propose the following 
metric tensor (( , }} for differential forms, 

(((a^^v; fc )H^I ( tf )(/? ^ fc | fc ^ 

W I W I 

where the \V a k \ is the factor arising from integrating the volume-form over V a k , and 

(a k ,a k ) (f3 k ,(T k ) 
\a k \ \a k \ 

is what we would expect for ((a k ,f3 k )), if the forms a k and f3 k were constant on a k , which is the product of 
the average values of a k , and (3 k . 

If we adopt this as our definition of the metric tensor for forms, we can recover the definition we obtained 
in the Hodge star operator. Starting from the definition from the smooth theory, 

({a k ,p k ))v = J a k /\*[3 k , 

and expanding this in terms of the metric tensor for discrete forms, and the primal-dual wedge operator, we 
obtain 

(((a k ,(3 k ))v,V ak ) = (a k A*[3 k ,V ak ), 

{a k ,n k ){(3 k ,o- k ) \V a u\ k k k k 
\a K \ \a K \ \a K \\*a K \ 
When we eliminate common factors from both sides, we obtain the expression, 



1 ,„*. ,, 1 



which is the expression we previously obtained in Definition \6.1\ of |J3 



The L 2 norm for discrete differential forms is given below. 



Definition 12.3. Given a primal discrete k-form a k e Q^K), its discrete L 2 norm is given by 

\\a k f d = {a k ,a k ) d 

a k £K 

1 ^ I * gfc | /„k k\2 

n \a k \ 

a k GK 1 1 
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Giv en these definitions, we can now reproduce some computations that were originally shown in lCastrill6n-L6ne3 
|2003j . 

Harmonic Functions. Harmonic functions <f> : M — > R can be characterized in a variational fashion as 
extremals of the following action functional, 

S^) = \f ||d^|| 2 v, 

where v is a Riemannian volume- form in M. The corresponding Euler-Lagrange equation is given by 

*d * d<j) = —A(f> = 0, 

which is the familiar characterization of harmonic functions in terms of the Laplace-Beltrami operator. 

The discrete action functional can be expressed in terms of the L 2 norm we introduced above for discrete 
forms, 



1 v-^ I* " 1 1 , 1 , - 

£ E 



2n \a 

The basic variations needed for the determination of the discrete Euler-Lagrange operator are obtained from 
variations that vary the value of the function <p at a given vertex vq, leaving the other values hxed. These 
variations have the form, 

(f) e = <f> + erj, 

where f) G QP(M; R) is such that (fj, vq) = 1, and (fj, v) = 0, for any v G — {vq}- This family of variations 
is enough to establish the variational principle. That is, we have 

S d (<f>s) 

E=0 



de 



n * — ' \<7 



(12.1) 



1 \~ k(jl \ 

- V Vii W.ir^sgQ ij 1 ;^ 



where sgn(cr 1 ;w) stands for the sign of a 1 with respect to v. Which is to say, sgn(cr 1 ;w) = 1 if a 1 — [v',v], 
and sgn(cr 1 ;u) = —1 if a 1 = [v,v']. On the other hand, 

(*d * d<j>, Vq) = | 1 | (d * d(p,*v ) 
(*d<t),d*v a ) 



\*v 
1 



— - V (*d0,*cr 1 }sgn(cr 1 ;uo) 
-Tl E ^(d^sgn^o), 

*V()\ \<T L \ 



•wo, 



where in the second to last equality, one has to note that the border of the dual cell of a vertex vq consists, 
up to orientation, in the dual of all the 1-simplices starting from vq. This is illustrated in Figure [2~T1 and 
follows from a general expression for the boundary of a dual cell that was given in Definition 15. 81 
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(a) Vertex v (b) *v (c) d * v 




(d) cr 1 y v (e) Vrcr 1 

Figure 21. Boundary of a dual cell. 

The sign factor comes from the relation between the orientation of the dual of the 1-simplices and that of 
d * vq. From this, we conclude that the variational discrete equation, given in Eauation ll2.ll is equivalent 
to the vanishing of the discrete Laplace-Bcltrami operator, 

*d * dcj) = -A = 0. 



Maxwell Equations. We can formulate the Maxwell equations of electromagnetism in a covariant fashion 
by considering the 1-form A (the potential) as our fundamental variable in a Lorentzian manifold X. The 
action functional for a Lagrangian formulation of electromagnetism is given by, 

S(A) = 1 f lld^jV 
•J x 

where || • || is the norm on forms induced by the Lorentzian metric on X, and v is the pseudo-Riemannian 
volume-form. The 1-form A is relate d to the 4-vecto r potential encountered in the relativistic formulation 
of electromagnetism (see, for example I Jackson! 1998j ). 



The Euler-Lagrange equation corresponding to this action functional is given by 

*d * dA = . 

In terms of the field strength, F — dA, the last equation is usually rewritten as 

dF = , *d*F = 0, 

which is the geometric formulation of the Maxwell equations. 

For the purposes of simplicity of exposition, we consider the special case where the Lorentzian manifold 
decomposes into X = M x R, where (M, g) is a compact Riemannian 3-manifold. In formulating the discrete 



12 



MATHIEU DESBRUN, ANIL N. HIRANI, MELVIN LEOK, AND JERROLD E. MARSDEN 



version of this variational problem, we need to generalize the notion of a discrete Hodge dual to take into 
account the pseudo-Riemannian metric structure. This can be subtle in practice, and to overcome this, we 
consider a special family of complexes instead. 

Let K' be a simplicial complex modelling M. For the sake of simplicity we consider M = R 3 although 
this is not strictly necessary. We now consider a discretization {i„}„gz of R. We define the complex K, 
modelling X = R 4 , the cells of which are the sets a — a' x {t n } C R 3 X R, and a — a' x (t n ,t n+ i) C R 3 x R 
for any a' G K' and n G Z. Of course, this is not a simplicial complex but rather a "prismal" complex, as 
shown in Figure 1221 




Figure 22. Prismal cell complex decomposition of space-time. 

The advantage of these cell complexes is the existence of the Voronoi dual. More precisely given any 
prismal cells a' x {t n } G K and a' x (t n ,t n+ i) G K, the Lorentz orthonormal to any of its edges coincide 
with the Euclidean one in R 4 and the existence of the circumcenter is thus guaranteed. In other words, the 
Lorentz circumcentric dual *K to K is the same as the Euclidean one in R 4 . 

Remark 12.2. Much of the construction above can be carried out more generally by considering arbitrary 
cell complexes in R 4 that are not necessarily prismal, as long as none of its 1-cells are lightlike. This causality 
condition is necessary to ensure that the circumcentric dual complex is well-behaved. However, it is sufficient 
for computational purposes that the complex is well- centered, in the sense that the Lorentzian circumcenter 
of each cell is contained inside the cell. These issues will be addressed in future work. 

Recall that the Hodge star * is uniquely defined by satisfying the following expression, 

a A */3 = {{a, /3))v , 

for all a, (3 G Vl k (X). The upshot of this is that the Hodge star operator depends on the metric, and since 
we have a pseudo-Riemannian metric, there is a sign that is introduced in our expression for the discrete 
Hodge star (Definition Itj.lfl that depends on whether the cell it is applied to is either spacelike or timelike. 
The discrete Hodge star for prismal complexes in Lorentzian space is given below. 

Definition 12.4. The discrete Hodge star for prismal complexes in Lorentzian space is a map 

* : VL^{K) — > il^(*K) defined by giving its action on cells in a prismal complex as follows, 

where \ ■ \ stands for the volume and the causality sign n{a k ) is defined to be +1 if all the edges of a k are 
spacelike, and —1 otherwise. 

The causality sign of 2-cells in a (2 + l)-space-time is summarized in Table 0] We should note that the 
causality sign for a 0-simplex, is always 1. This is because a 0-simplex has no edges, and as such the 

statement that all of its edges are spacelike is trivially true. 
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W 1 ) 








-1 



-1 



-1 



This causality term in the discrete Hodge star has consequences for the expression for the discrete norm 
(Definition 112.31) . which is now given. 

Definition 12.5. Given a primal discrete k-form a k S il^(K), its discrete L 2 Lorentzian norm is given 
by, 



a 



k i|2 



I Lor.d 



= - V (a k ,a k )(*a k ,*a k ) 
a k ei< 



a k eK 



Having defined the discrete Lorentzian norm, we can express the discrete action as 



S d (A) 



1, 



Idyll 



Lor.d 



i J2 (dA,a*){*dA,*a 2 ) 

a 2 eK 

5 E "(OKW,^ 



o 2 eK 1 1 

The basic variations needed to determine the discrete Eulcr-Lagrange operator are obtained from variations 
that vary the value of the 1-form A at a given 1-simplex <7q, leaving the other values fixed. These variations 
have the form, 

A £ = A e + efj, 

where fj € Cl\(K) is given by (fj, Oq) = 1 for a fixed interior <Tq S A' and (fj, a 1 ) = for a 1 ^ ctq. The 
derivation of the variational principle gives 



d_ 

ds 



S M = \ E if^^XdA^Kd^ 2 



1 E 



★ (7 



U(<r 2 )(dA, a 2 )(dr),a 2 ) 



■K(o 2 )(dyl, a 2 )sgn(a 2 ,a,i), 
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which vanishes for all the basic variations above. On the other hand, we now expand the discrete 1-form 
*d * dA. For any <Jq £ K, we have that 



*d*dA,o%) = J-^yK(<rJ) (d*d A, ★<t 1 } 
:K((rJ){*di, d-kal) 



jfii 
jfii 



-«(°o) H sgn((i 2 ,cri)(*dA,* ( j 2 ) 



where the sign sgn(er 2 , a 1 ) stands for the relative orientation between a 2 andc 1 . Which is to say, sgn(tr 2 , a 1 ) — 
1 if the orientation induced by a 2 on a 1 coincides with the orientation of cr 1 , and sgn(<T 2 ,cr 1 ) = —1 other- 
wise. For the second to last equality, one has to note that the border of the dual cell of an edge Oq consists, 
conveniently oriented with the sgn operator, of the union of the duals of all the 2-simplices containing <Jq. 
This statement is the content of Definition 15.81 which gives the expression for the boundary of a dual cell, 
and was illustrated in Figure EHI for the case of n-dimensional dual cells. 

By comparing the two computations, we find that for an arbitrary choice of (Tq G K, (*d * dA, <7q) is 
equal (up to a non-zero constant) to 5Sd(A), which always vanishes. It follows that the variational discrete 
equations obtained above is equivalent to the discrete Maxwell equations, 

*d * dA = 0. 



13. Extensions to Dynamic Problems 

It is desirable to leverage the exactness properties of the operators of discrete exterior calculus to construct 
numerical algorithms with discrete conservation properties. For these purposes, it is appropriate to extend 
the scope of DEC to incorporate dynamical behavior, by addressing the issue of discrete diffeomorphisms 
and flows. 

As discussed in the previous section, DEC and discrete mechanics have interesting synergistic proper- 
ties, and in this section we will explore a groupoid interpretation of discrete mechanics that is particularly 
appropriate to formulating the notion of pull-back and push-forward of discrete differential forms. 



13.1. Groupoid Interpretation of Discrete Variational Mechanics. The groupoid formulation of dis- 
crete mechanics is particularly fruitful and natural, and it serves as a unifying tool for understanding the 
variational for mulation of discre te L agrangian mechanics, and di screte Euler-Poincare reduction, as discussed 
in the work of IWeinsteinl |l99fij and lMarsden et all |l999ll200Cj . 

The groupoid interpretation of discrete mechanics is most clearly illustrated if we consider the discretiza- 
tion of trajectories on TQ in two stages. Given a curve 7 : R + — > TQ, we consider a discrete sampling given 

by 

9l = y(ih) e TQ. 

We then approximate TQ by Q x Q, and associate to gi two elements in Q. We denote this by 

9i ^ {q°,ql)- 
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Or equivalently, in the language of groupoids, see ICannas da Silva and Weinsteinl 19991 ]. IWeinsteinl 2001], 
we have 

G 

■■ \ 1 ■ 



where a is the source map, and f3 is the target map. Then, 

9i i-> (a(gi),/3(gi)) = {q°,ql) 

This can be visualized as 




ii = a i9i) q} = Pirn) 

A product • : G^ — > G is defined on the set of composable pairs, 

G {2) :={(g,h)£GxG\{3(g) = a(h)}. 

The groupoid composition g ■ h is defined by 

a(g ■ h) = a(g), 
0{g ■ h) = (3(h). 

This can be represented graphically as follows, 



g-h 



a(g) = a(g ■ h) /3(g) = a(h) (3(h) = f3(g ■ h) 

The set of composable pairs is the discrete analogue of the set of second-order curves on TQ. A curve 
7 : WL + — > TQ is said to be second-order if there exists a curve q : R + — > Q, such that, 

7(*) = (?(*), «(*))• 

The corresponding condition for discrete curves is that given a sequence of points in QxQ, (gj, ■ ■ ■ , (<?p, 9p), 
we require that 

Qi = Q°+v 

This implies that the discrete curve on Q x Q is derived from a (p+ l)-pointed curve (goi • • • j 9p) on Q, where 

?i= /«i+i' if 0<*<p; 

iii=p- 

This condition has a direct equivalent in groupoids, 

P(9i) = Qi = Qi+i = a{g l+ i). 
Which is to say that the sequence of points in Q x Q are composable. In general, this hierarchy of sets is 
denoted by 

G (ri :={(9i,...,g P )€Gi> \f3{ gi ) = a{g i+l )} : 

where G^ ~ Q. 
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air 1 ) =/?(<?) 



a 



Visualizing Groupoids. In summary, composition of groupoid elements, and the inverse of groupoid 
elements can be illustrated by Figure 1231 As we will see in the next subsection, representing discrete 



diffeomorphisms as pair groupoids is the natural method of ensuring that the mesh remains nondegenerate. 

13.2. Discrete Diffeomorphisms and Discrete Flows. We will ado pt the point of view of rep resenting 
a discrete diffeomorphism as a groupoid, which was first introduced in IPekarskv and Westl [2003J, and ap- 
propriately modify it to reflect the simplicial nature of our mesh. In addition, we will address the induced 
action of a discrete diffeomorphism on the dual mesh. 

Definition 13.1. Given a complex K embedded in V , and its corresponding abstract simplicial complex M , 
a discrete diffeomorphism, ip 6 Diff f j(M), is a pair of simplicial complexes K\, K^, which are realizations 
of M in the ambient space V. This is denoted by tp(M) = (Kx^Kz). 

Definition 13.2. A one-parameter family of discrete diffeomorphisms is a map ip : I — > Diff^M), 




Figure 23. Groupoid composition and inverses. 



such that, 



Vs,i G I. 
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Since we are concerned with evolving equations represented by these discrete diffeomorphisms, and mesh 
degeneracy causes the numerics to fail, we introduce the notion of non-degenerate discrete diffeomorphisms, 

Definition 13.3. A non-degenerate discrete diffeomorphism ip = {K\,K%) is such that K\ and K 2 
are non- degenerate realizations of the abstract simplicial complex M in the ambient space V . 

Notice that it is sufficient to define the discrete diffeomorphism on the vertices of the abstract complex 
M, since we can extend it to the entire complex by the relation 

<p{[V0, -,Vk\) = ([7Tl^(uo), ...,-Klip(Vk)}, [ir 2 (p(vo), ...,W 2 <p(vk)]). 

If X £ is a material vertex of the manifold, corresponding to the abstract vertex w, that is to say, 

TTi<Pt(w) — X,\ft £ /, the corresponding trajectory followed by X in space is x — n 2 ip t (w). Then, the 
material velocity V(X,t) is given by 

dn 2 tp s {w) 



V(tti («>),*) 



ds 

and the spatial velocity v(x,t) is given by 

dL Ps (tp- 1 ( x )) 



v{n 2 (w),t) = V(m(w),t) 



ds 



s=t 

The distinction between the spatial and material representation is illustrated in Figure [ 




The material velocity field can be thought of as a discrete vector field with the vectors based at the vertices 
of K, which is to say that Tip t £ Xd(K), is a discrete primal vector field. Notice that ipt on K induces a 
map -k(pt on the vertices of the dual *K , by the following, 

■kcp t (c[v , . . .,v n ]) = {c[ni(p f (v ), ■ ■ ■ ,7ri<pt(u n )],c[7r2^t(«o), ■ ■ ■ , ^2ft(v n )])- 
Similarly then, T-ktpt £ Xd(*K) is a discrete dual vector field. 

Comparison with Interpolatory Methods. At first glance, the groupoid formulation seems like a 
cumbersome way to define a one-parameter family of discrete diffeomorphisms, and one may be tempted to 
think of extending ip t to the ambient space. We would then be thinking of ipt : V — > V. This is undesirable 
since given cpt and ip a which are non-degenerate flows, their composition tp t o^ S) which is defined, may result 
in a degenerate mesh when applied to K . Thus, non-degenerate flows are not closed under this notion of 
composition. 
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If we adopt groupoid composition instead at the level of vertices, we can always be sure that if we compose 
two nondegenerate discrete diffeomorphisms, they will remain a nondcgenerate discrete diffcomorphism. 

Discrete Diffeomorphisms as Pair Groupoids. The space of discrete diffeomorphisms naturally has 
the structure of a pair groupoid. The discrete analogue of TDiff(M) from the point of view of temporal 
discretization is the pair groupoid Diff(M) x Diff(M). In addition, we discretize Diff(M) using Diff^M), 
which is in turn a pair groupoid involving realizations of an abstract simplicial complex in an ambient space. 

13.3. Push-Forward and Pull-Back of Discrete Vector Fields and Discrete Forms. For us to 

construct a discrete theory of exterior calculus that admits dynamic problems, it is critical that we introduce 
the notion of push-forward and pull-back of discrete vector fields and discrete forms under a discrete flow. 

Push-Forward and Pull-Back of Discrete Vector Fields. 

The push-forward of a discrete vector field satisfies the following commutative diagram, 

x ■ -N 

Tf 
N 

and the pull-back satisfies the following commutative diagram, 

K 
f 



By appropriately following the diagram around its boundary, we obtain the following expressions for the 
push-forward and pull-back of a discrete vector field. 

Definition 13.4. The push-forward of a dual discrete vector field X G Xd{*K), under the map 
f : K —> L, is given by its evaluation on a dual vertex <r = *er™ G 

f.X(*v n )=Tf.X{*{f-\v n ))). 

Definition 13.5. The pull-back of a dual discrete vector field X G Xd(*L), under the map f : K — > L, 

is given by its evaluation on a dual vertex &q = •ka n G (*K)^ , 

f*X{*o- n ) = (f-%X(*a n ) = TCr 1 ) • X(*(f(a n ))). 

Pull-Back and Push-Forward of Discrete Forms. A natural operation involving exterior calculus in 
the context of dynamic problems is the pull-back of a differential form by a flow. We define the pull-back of 
a discrete form as follows. 

Definition 13.6. The pull-back of a discrete form a k G il^(L), under the map f : K — > L, is defined 
so that the change of variables formula holds, 

(ra k ,a- k ) = (a k ,f(o- k )), 

where a k G K . 

We can define the push- forward of a discrete form as its pull-back under the inverse map as follows. 

Definition 13.7. The push-forward of a discrete form a k G Sl^K), under the map f : K — > L is 

defined by its action on a k G L, 

(f*a k ,a k ) = <(/- 1 )*a fc ,a fc ) = (a k , f-\a k )). 
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Naturality under Pull-Back of Wedge Product. We find that the discrete wedge product we introduced 
in [JE1 is not natural under pull-back, which is to say that the relation 

r(aA(3) = f*aAf*l3, 

does not h old in general. Howeve r, a metric independent definition that is natural under pull-back was 
proposed in ICastrillon-Lopezl |2003| . 



Definition 13.8 l|Castrill6n- Given a primal discrete k-form a k G il^(K), and a primal 

discrete I- form (3 l G fl l d (K), the natural discrete primal-primal wedge product, A : il d (K) x il l d (K) — » 
fl d +l (K), is defined by its evaluation on a (k + l)-simplex a k+l — [vq, . . . ,Vk+i], 

(o; fc A/3V fc +')= * + ]T sign(r)a^P(r(a k+l )). 

In contrasting this definition to that given by Definition 18.11 we see that the geometric factor 

\a k+l n*v T(k) \ 
\a k + l \ 

has been replaced by 



k + l + 1 

in this alternative definition. By replacing the geometric factor which is metric dependent with a constant 
factor, Definition 1 1 3 . 81 becomes natural under pull-back. 

The proofs in [JHlthat the discrete wedge product is anti-commutative, and satisfies a Leibniz rule, remain 
valid for this alternative discrete wedge product, with only trivial modifications. As for the proof of the 
associativity of the wedge product for closed forms, we note the following identity, 

V k fc+ ' n ^r( fc )l _ V 1 _ (k | n, 

W k+l \ ~ k + i + i + 

r£S k+l + 1 1 1 TES k+l + 1 

which is a crucial observation for the original proof to apply to the alternative wedge product. 

14. Remeshing Cochains and Multigrid Extensions 

It is sometimes desirable, particularly in the context of multigrid, multiscale, and multiresolution compu- 
tations, to be able to represent a discrete differential form which is given as a cochain on a prescribed mesh, 
as one which is supported on a new mesh. Given a differential form uj k G il k (K), and a new mesh M such 
that \K\ — \M\, we can define it at the level of cosimplices, 

Vr fe GMW, (u k ,r k }= J2 sgn(r fc ,a fc ) lK ; T/ n ^' (^,a fc ), 

er fc eif( fc ) ' a 

and extend this by linearity to cochains. Here, sgn(r fe , a k ) is +1 if the orientation of r k and a k are consistent, 
and —1 otherwise. Since fc-skeletons of meshes that are not related by subdivision may not have nontrivial 
intersections, intersections of support volumes are used in the remeshing formula, as opposed to intersections 
of the fc-simplices. 

We denote this transformation at the level of cochains as, Tk,m ■ C k (K) — » C k (M). This has the natural 
property that if we have a fc-volume U that can be represented as a chain in either the complex K or the 
complex M, that is to say, U k = a k + . . . + a k = r k + . . . + r k , then we have 

m m \v k nv fc 

(^,^ + ... + r^) = ^(c,rf>=^ sgn(T*> fc ) ' T - J \ u, k ,a k ) 
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m 1 W k l~l V k\ 

i=l 3=1 1 °j 1 

3=1 

Which is to say that the integral of the differential form over U k is well-defined, and independent of the 
representation of the differential form. 

Note that, in particular, if we choose to coarsen the mesh, the value the form takes on a cell in the coarser 
mesh is simply the sum of the values the form takes on the old cells of the fine mesh which make up the new 
cell in the coarser mesh. 

Non-Flat Manifolds. The case of non-flat manifolds presents a challenge in remeshing akin to that 
encountered in the discretization of differential forms. In particular, if the two meshes represent different 
discretizations of a non-flat manifold, they will in general correspond to different polyhedral regions in the 
embedding space, and not have the same support region. 

We assume that our discretization of the manifold is sufficiently fine that for every simplex, all its vertices 
are contained in some chart. Then, by using these local charts, we can identify support volumes in the 
computational domain with n-volumes in the manifold, and thereby make sense of the remeshing formula. 



15. Conclusions and Future Work 

We have presented a framework for discrete exterior calculus using the cochain representation of discrete 
differential forms, and introduced combinatorial representations of discrete analogues of differential operators 
on discrete forms and discrete vector fields. The role of primal and dual cell complexes in the theory are 
developed in detail. In addition, extensions to dynamic problems and multi-resolution computations are 
discussed. 

In the next few paragraphs, we will describe some of the future directions that emanate from the current 
work on discrete exterior calculus. 

Relation to Computational Algebraic Topology Since we have introduced a discrete Laplace-deRham 
operator, one can hope to develop a discrete Hodge-deRham theory, and relate the deRham cohomology of 
a simplicial complex to its simplicial cohomology. 

Extensions to Non-Flat Manifolds. The intrinsic notion of what constitutes the discrete tangent space 
to a node on a non-flat mesh remains an open question. It is possible that this notion is related to a choice 
of discrete connection on the mesh, and it is an issue that deserves further exploration. 

Generalization to Arbitrary Tensors. The discretization of differential forms as cochains is particularly 
natural, due to the pairing between forms and volumes by integration. When attempting to discretize 
an arbitrary tensor, the natural discrete analogue is unclear. In particular, while it is possible to expand 
an arbitrary tensor using the tensor p roduct of covaria nt and contravariant one-tensors, this would be 
cumbersome to represent on a mesh. In lLeok et al.l 2fJ03j, which is on discrete connections, we will see Lie 
group-valued discrete 1-forms, and one possible method of discretizing a (p, g)-tensor that is alternating in 
the contravariant indices, is to consider it as a (0, g)-tensor-valued discrete p-form. 
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It wo uld be particu larly interesting to explore this in the context of the elasticity complex (see, for 



it wo uld De particuia 
example, lArnoldl [2002J ) , 



se(3)C yC°°(n,R 3 ) —^C°°(Q,S) — T -^C°°(Q,S) ^^C 00 ^,^ 3 ) >0, 

where § is the space of 3 x 3 symmetric matrices. One approach to discreti ze this was suggested in lArnoldl 
2002J, which cites the use of the Bernstein-Gelfand-Gelfand resolution in lEastwoodl |2000j to derive the 
elasticity complex from the deRham complex. Alternatively, it might be appropriate in the context of the 
elasticity complex to consider Lie algebra-valued discrete differential forms. 

Convergence and Higher-Order Theories. The natural question from the point of view of numerical 
analysis would be to carefully analyze the convergence properties of these discrete differential geometric 
operators. In addition, higher-order analogues of the discrete theory of exterior calculus are desirable from 
the point of view of computational efficiency, but the cochain representation is attractive due to its conceptual 
simplicity and the elegance of representing discrete operators as combinatorial operations on the mesh. 

It would therefore be desirable to reconcile the two, by ensuring that high-order interpolation and com- 
binatorial operations are consistent. As a low-order example, Whitney forms, which are used to interpolate 
differential forms on a simplicial mesh, have the nice property that taking the Whitney form associated with 
the coboundary of a simplicial cochain is equal to taking the exterior derivative of the Whitney form asso- 
ciated with the simplicial cochain. As such, the coboundary operation, which is a combinatorial operation 
akin to finite differences, is an exact discretization of the exterior derivative, when applied to the degrees of 
freedom associated to the finite-dimensional function space of Whitney forms. 

It would be interesting to apply subdivision surface techniques to construct interpolatory spaces that 
are compatible with differential geometric operations that are combinatorial operations on the degrees of 
freedom. This will result in a massively simplified approach to higher-order theories of discrete exterior 
calculus, by avoiding the use of symbolic computation, which would otherwise be necessary to compute the 
action of continuous exterior differential operators on the polynomial expansions for differential forms. 
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